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ENVIRONMENTAL  INFLUENCES  IN  THE  SIMULATION 
OF  A  SOLAR  SPACE  HEATING  SYSTEM 

CHAPTER  I 

INTRODUCTION 

Background 

After  about  35  years  of  limited  development,  solar  energy 
has  recently  become  a  viable  and  economically  attractive  energy  sup¬ 
plement,  capable  of  reducing  our  reliance  on  increasingly  scarce  fos¬ 
sil  fuel.  Complete  systems  are  now  being  sold  and  installed  that 
have  a  storage  medium  with  the  capability  to  provide  from  one  to  three 
days  space  heating  requirements. 

Such  systems  involve  a  mix  of  both  solar  and  conventional 
heating  systems,  since  a  solar  system  which  provides  all  of  the  build¬ 
ing  heat  demands  is  normally  neither  economically  nor  architecturally 
attractive.  Given  the  short  term  storage  and  the  unpredictable  nature 
of  the  solar  input,  the  size  of  the  conventional  system  in  this  com¬ 
bination  must  be  sufficient  to  provide  all  of  the  energy  demands  of 
the  building  under  the  most  adverse  conditions  expected.  The  size  of 
the  solar  system,  on  the  other  hand,  is  determined  primarily  by  the 
marginal  cost  of  the  collector/storage  system  in  comparison  with  the 


cost  of  the  conventional  fuel  saved. 
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The  amount  of  fuel  saved  is,  in  turn,  a  function  of  the  per¬ 
formance  characteristics  of  the  solar  system  and  the  heat  loss  charac¬ 
teristics  of  the  building.  In  many  studies  published  in  recent  years, 
the  optimum  collector  size,  in  an  economic  sense,  has  been  determined 
from  numerical  simulation  of  hourly  collector  performance  and  building 
heat  demand  over  a  representative  period.  The  usual  method  of  model¬ 
ing  building  heat  loss  in  these  studies  has  taken  into  account  temper¬ 
ature  effects  only.  A  value  for  the  building  heat  loss  per  hour  per 
unit  temperature  difference  is  determined  using  the  procedures  out¬ 
lined  in  the  ASHRAE  Handbook  of  Fundamentals  [2]  for  the  building  con¬ 
struction  characteristics  and  climatological  design  parameters  at  the 
location  under  study.  It  is  this  value  and  the  hourly  values  of  air 
temperature  which  then  determine  the  building  heat  loss  in  the  simula¬ 
tion.  In  actuality,  these  hourly  values  of  building  heat  loss  will 
vary  not  only  with  changes  in  air  temperature,  but  with  the  wind 
through  its  effects  on  convective  losses  and  air  infiltration  into 
the  structure,  and  with  solar  radiation  effects  on  the  exterior  of 
the  structure. 


Qbj  ective 

Modeling  building  heat  loss  as  a  function  of  temperature 
only,  in  effect  the  degree  day  method,  has  been  shown  [20]  to  be  rea¬ 
sonably  accurate  over  a  long  period  of  time,  but  does  not  necessarily 
accurately  simulate  the  short  term  nature  of  building  heat  loss.  The 


objective  of  this  study  is  to  determine  whether  these  hourly  depar¬ 
tures  in  building  heat  loss  from  the  values  given  by  a  simplified 
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degree  day  approach  have  a  significant  effect  on  the  selection  of  an 
optimum  collector  size  for  the  solar  system.  In  order  to  pursue  the 
objective,  a  numerical  model  was  developed  which  simultaneously  sim¬ 
ulates  the  hourly  interactions  between  the  solar  system,  the  house, 
and  the  surrounding  environment,  using  both  the  simplified  model  and 
an  expanded  parameter  model  which  includes  the  effects  of  air  movement 
over  the  structure  and  the  influence  of  radiation  on  the  structure  it¬ 
self.  Using  the  results  from  this  simulation  an  optimum  collector 
size  was  determined  from  the  energy  requirements  given  by  each  model 
and  a  comparison  made  between  the  two  approaches. 
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CHAPTER  II 

SOLAR  RADIATION  AND  METEOROLOGICAL  DATA 

The  model  developed  for  this  study  is  designed  to  use  a 
SOLMET  [24]  data  tape  for  the  input  values  of  solar  radiation  and 
meteorological  parameters.  These  data  tapes,  obtained  from  the 
National  Climatic  Center,  contain  hourly  solar  radiation  data  and 
meteorological  observations. 

The  radiation  values  on  these  tapes  are  of  two  types.  At 
stations  with  solar  observation  facilities,  they  are  the  observed 
values  corrected  for  known  instrumentation  errors  [24].  For  those 
stations  without  solar  observations,  radiation  values  were  modeled 
through  a  regression  procedure  based  on  local  meteorological  con¬ 
ditions  and  radiation  values  at  the  nearest  station  with  observed  val 
ues  [24].  These  regression  modeled  values  were  originally  reported 
to  be  accurate  within  a  standard  error  of  +/-  1.7%  per  month  when 
compared  with  observed  values  [24].  However,  these  values  have  re¬ 
cently  been  questioned  by  Hoyt  [11]  who  reports  the  monthly  errors 
may  be  as  large  as  7%  at  some  stations.  The  SOLMET  values  are  still 
apparently  the  most  accurate  estimate  of  hourly  solar  radiation,  par¬ 
ticularly  when  compared  with  the  frequently  used  ASHRAE  clear  sky  val 
ues  which  Felske  [7]  reports  give  large  overestimates  of  available 
solar  energy. 
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The  SOLMET  tapes  with  regression  modeled  solar  data  have 
values  only  for  global  radiation  on  a  horizontal  surface.  These  val¬ 
ues  must  be  separated  into  direct  and  diffuse  components  for  solar 
system  performance  calculations  and  this  separation  was  accomplished 
using  the  method  developed  by  Liu  and  Jordan  [14]  and  modified  by 
Hunn  [12]  for  computer  analysis.  In  this  method,  Liu  and  Jordan's 
attenuation  index,  KT,  is  taken  to  be 

KT=A*(1  -  . 1*CC) 

where  A  is  the  transmissivity  and  CC  the  total  opaque  cloud  cover  in 
tenths.  The  diffuse  component  is  separated  from  the  global  horizon¬ 
tal  value  through  a  correlation  of  Liu  and  Jordan's  value  for  the  ra¬ 
tio  of  diffuse  to  horizontal  radiation,  D/H,  with  the  attenuation  in¬ 
dex,  KT,  by  assuming  this  ratio  can  be  applied  to  the  hourly  SOLMET 
values : 

(D/H) =1-1. 3575 *KT 
RHDU= (D/H) *RHT 

where  RHDU  is  the  diffuse  component  and  RHT  the  global  horizontal  ra¬ 
diation.  The  horizontal  direct  value  is  then  obtained  by  subtracting 
the  diffuse  from  the  global  and  a  direct  normal  value  calculated  from 
this  value  and  the  solar  altitude: 

RD=RHT-RHDU 

RDN=RD/sin(ct) 

where  RD  is  the  horizontal  direct  value,  RDN  the  direct  normal  value, 
and  alpha  the  solar  altitude. 

Those  SOLMET  stations  which  have  corrected  observed  values 
also  have  values  of  direct  normal  radiation  on  the  tape  [24].  The 
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computer  model  developed  for  this  study  is  designed  to  take  either 
type  tape  as  input,  and,  depending  upon  an  input  flag,  use  either  the 
given  direct  normal  value  or  one  calculated  as  indicated  above.  The 
other  parameters  from  the  SOLMET  tapes  used  in  this  study  are  the 
year,  month,  day,  solar  time,  temperature,  wind  speed  and  direction, 
and  total  opaque  cloud  cover. 

The  SOLMET  data  tape  for  Oklahoma  City,  with  13  years  (1952- 
1964)  of  hourly  radiation  and  hourly  meteorological  data  was  used  for 
this  study.  To  reduce  the  computer  time  necessary  for  the  simulation, 
the  13  year  data  set  was  used  to  produce  an  average  year  which  consist¬ 
ed  of  the  13  year  average  hourly  value  for  each  parameter  for  each  of 
the  8760  hours  in  a  365  day  year.  The  average  year  values  were  trans¬ 
ferred  to  disk  storage  for  more  ready  accessibility  and  served  as  the 
data  set  for  all  computer  simulations.  Beckman,  et  al.  [3]  have  shown 
that  the  use  of  an  average  year  in  the  simulation  of  solar  system 
operation  closely  approximates  the  results  obtained  by  using  the  full 
data  set  and  averaging  the  results,  thus  providing  an  accurate  esti¬ 
mate  of  long  term  performance. 
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CHAPTER  III 

HOUSE  MODELS 

Construction  Characteristics 

The  houses  considered  in  this  study  were  modeled  as  single 
story  residences  with  an  unventilated  attic  beneath  a  peaked  roof, 
and  a  full  basement,  as  illustrated  in  Figure  1,  in  which  the  solar 
system  storage  vessel  and  conventional  heating  system  were  placed. 

The  floor  space  for  the  three  different  house  sizes  considered,  de¬ 
noted  as  size  I,  II,  and  III,  is  93  square  meters,  135  square  meters, 
and  180  square  meters.  Each  house  faces  south  and  is  taken  to  be  a 
rectangular  structure  with  width  to  length  ratio  of  about  1:1.5  as 
suggested  by  Olgay  [21]  for  optimum  winter  thermal  performance  in 
temperate  to  cool  climates.  The  solar  collector  array  is  taken  as  an 
integral  part  of  the  south  facing  roof  and  has  a  slope  from  the  hori¬ 
zontal  of  latitude  plus  15  degrees,  the  value  reported  by  Liu  and 
Jordan  [14]  to  be  optimum  for  solar  space  heating.  Since  the  model 
is  designed  to  run  for  different  locations  and  many  collector  areas, 
the  remaining  dimensions  of  the  roof,  such  as  the  slope  and  area  of 
the  north  roof,  the  area  of  the  attic  walls,  and  the  height  of  the 
peak,  are  calculated  in  the  program  when  the  collector  tilt  and  area 
of  collector  to  be  integrated  into  the  roof  are  specified.  Each  house 
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Illustration  of  house  model. 


has  a  total  window  area  of  approximately  10%  of  the  house  wall  area, 
and  following  Olgay  [21],  the  distribution  is  weighted  in  favor  of 
southern  exposure  for  maximum  solar  gain  in  winter.  The  window  area 
distribution  is  set  at  60%  of  total  window  area  on  the  south  wall, 

10%  on  the  east  and  west  walls,  and  20%  on  the  north  wall.  Each  house 
has  two  entrances,  one  on  the  north  side  and  one  on  the  south  side  of 
the  structure.  The  dimensions  for  each  house  size  are  given  in  Table 
A-l  of  Appendix  A. 

Insulation 

For  each  house  size  in  the  model,  three  different  insulation 
and  weatherstripping  levels  were  considered,  corresponding  to  poorly 
insulated,  moderately  insulated,  and  well  insulated  houses,  and  de¬ 
noted  as  types  I,  II,  and  III.  Type  1  consists  of  minimum  wall  and 
ceiling  insulation,  single  pane  windows  with  aluminum  sash,  no  storm 
windows  or  doors,  and  minimum  weatherstripping.  Type  II  has  moderate 
wall  and  ceiling  insulation,  double  pane  windows  with  wood  sash,  storm 
doors,  and  moderate  weatherstripping.  Type  III  has  heavy  wall  and 
ceiling  insulation,  double  pane  windows  with  wood  storm  sash,  insulated 
doors  and  storm  doors,  and  is  tightly  weatherstripped.  The  construc¬ 
tion  characteristics  and  heat  transmission  coefficients  (U-values) 
for  each  component  of  the  three  insulation  levels  are  given  for  each 
type  in  Tables  A-2  and  A-3  of  Appendix  A. 
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CHAPTER  IV 

BUILDING  HEAT  LOSS  MODELS 

Introduction 

Calculation  of  the  energy  requirements  for  the  heating  of  a 
building  may  be  carried  out  at  a  very  simple  level  to  provide  approx¬ 
imate  results  or  at  increasingly  more  complex  levels  for  more  refined 
determination  of  system  performance. 

At  the  lower  end  of  the  scale  is  the  heat  loss  per  unit  tem¬ 
perature  difference  method,  as  outlined  in  the  ASHRAE  Handbook  of 
Fundamentals  [2].  This  method  is  based  upon  the  assumption  of  steady 
state  heat  transfer  through  the  building  fabric  as  a  function  of  in¬ 
side  to  outside  temperature  difference.  In  this  simplified  model, 
denoted  in  this  study  as  Model  I,  an  overall  building  heat  loss  co¬ 
efficient  calculated  using  the  methods  outlined  in  the  ASHRAE  Handbook 
provides  the  basis  for  determination  of  energy  requirements  for  space 
heating.  This  simplified  method  is  the  one  most  frequently  used  in 
engineering  calculations  and  has  served  as  the  model  for  building  heat 
loss  in  many  studies  of  solar  energy  for  space  heating,  for  example 
the  widely  quoted  study  by  Lof  and  Tybout  [15]. 

At  the  upper  end  of  the  scale  is  the  numerical  simulation 
of  building  heat  loss  by  the  rigorous  method  of  simultaneous  solution 
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of  the  non-steady  state  balance  equations  for  each  hour  for  each  com¬ 
ponent  of  the  structure.  This  approach  is  presented  in  reference  20, 
which  outlines  the  National  Bureau  of  Standards'  computer  program  for 
heating  and  cooling  loads,  a  sophisticated  model  primarily  used  to  per¬ 
form  design  studies  of  heating  and  cooling  systems  for  large  struc¬ 
tures.  The  degree  of  complexity  of  this  approach  is,  however,  very 
time  consuming  and  a  model  at  an  intermediate  level  was  developed  for 
this  study.  Denoted  here  as  Model  II,  it  takes  into  account  wind  and 
solar  radiation  effects  on  the  building,  as  wel'  as  temperature  dif¬ 
ferences,  and  will  be  used  to  determine  if  these  influences  are  impor¬ 
tant  considerations  in  the  simulation  of  a  solar  space  heating  system. 

Model  I 

The  heat  loss  from  a  house  consists  of  the  energy  lost  by 
transmission  through  the  walls,  basement,  ceiling,  and  other  exposed 
surfaces  and  by  energy  carried  away  by  exfiltrating  air.  The  deter¬ 
mination  of  an  overall  building  heat  transmission  coefficient,  the 
heating  load  imposed  by  infiltrating  air,  and  the  combination  of  the 
two  effects  into  an  overall  coefficient  of  heat  loss  per  unit  temper¬ 
ature  difference  are  presented  below. 

Transmission  Heat  Loss 

Transmission  heat  loss  through  a  given  building  surface, 
under  the  assumption  of  steady  state  conditions,  may  be  determined 
using  the  procedures  outlined  in  the  ASHRAE  Handbook  of  Fundamentals 
[2].  Using  a  heat  transmission  coefficient  (U-value)  for  the  con¬ 
struction  characteristics  of  the  surface,  the  total  heat  transmission 
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through  the  surface  is 

Q  -  U*A*(ti  -  to)  (1) 

Q  a  heat  transfer  rate,  kJ/h 
ti,  to  =  temperature  at  inside/outside  surface,  C 

A  =  total  area  of  the  surface,  sq  m 
U  *  heat  transmission  coefficient,  kJ/(sq  m-h-K) 

The  heat  transmission  coefficient  can  be  expressed  as 
U  -  1/R 

R  =  thermal  resistance  of  the  ith  element  in  the 
composite  structure,  given  by 

R  =  x/k  for  homogeneous  materials,  where 

x  ■  thickness  of  the  element,  m 

k  »  thermal  conductivity  of  the  element,  kJ/(sq  m-h-C) 

R  ■  1/C  for  air  spaces,  where  C  is  the  thermal  conductance 
R  3  1/hi,  1/ho  at  the  inside  and  outside  surfaces,  and 
h  =  the  surface  film  conductance. 

The  surface  film  conductance,  which  accounts  for  convective  and  ra¬ 
diative  losses  at  the  surface,  is  determined  experimentally  for  a 
given  surface  and  can  be  expressed  as 

h  =  he  +  hr 

The  term  he  is  an  experimentally  determined  convective  heat  transfer 

coefficient  and  hr  is  a  radiative  transfer  coefficient  between  the 

surface  and  its  surroundings,  given  by 

hr  -  a*Fa*Fb*(Tl4  -  T24)/(Tl  -  T2) 

o  *  Stefan  Boltzman  constant 

Fa  *  dimensionless  factor  to  account  for  the  shape  and 
orientation  of  the  surface  to  its  surroundings 
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Fb  =  dimensionless  factor  which  accounts  for  the  emission 
and  absorption  characteristics  of  the  surface  and  its 
surroundings 

Tl,  T2  =  temperature  of  surface  and  surroundings,  K 
The  radiation  coefficient,  hr,  is  about  1/5  the  magnitude  of  the  con¬ 
vective  component  of  the  surface  film  conductance  and  varies  little 
[6]  over  the  range  of  temperatures  of  interest,  having  only  a  2% 
change  with  a  10  degree  C  change  in  temperature.  A  constant  value  of 
14  kJ  per  sq  m  per  hr  per  degree  C  [2]  was  used  for  this  parameter  in 
this  study. 

For  a  surface  of  given  construction  and  thickness,  the  heat 
transmission  coefficient,  U,  may  be  calculated  for  steady  state  con¬ 
ditions  if  the  thermal  conductivities,  thicknesses,  and  surface  resis¬ 
tances  are  known.  The  ASHRAE  Handbook  contains  extensive  tables  of 
resistances,  conductivities,  and  transmission  coefficients  for  most 
building  construction  materials,  determined  from  a  vast  amount  of  ex¬ 
perimental  and  theoretical  work.  These  tables  were  used  to  determine 
the  U-values  for  the  building  components  considered  in  this  study. 

Using  the  U-values  determined  from  the  ASHRAE  tables,  the 
total  transmission  heat  loss  through  the  building  fabric  per  unit 
time  may  be  expressed  as 

Qt  =  (ti-to){UgAg  +  UdAd  +  UwAw  +  UcAc}  +  (ti-tg)UbP  (2) 
where  U  and  A  refer  to  the  heat  transmission  coefficients  and  surface 
areas,  P  is  the  perimeter  of  the  house,  and  the  subscripts  g,  d,  w,  c, 
b  refer  to  window  glass,  exterior  doors,  exterior  walls,  ceiling,  and 
basement,  respectively.  Ub  is  an  effective  heat  transmission  coeffi¬ 
cient  per  unit  length  for  the  basement  walls  and  floor,  determined  as 
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outlined  in  the  ASHRAE  Handbook  for  below  grade  basements,  and  the 
value  for  Uc  is  an  overall  coefficient  of  heat  transfer  for  the  com¬ 
bined  ceiling  and  roof.  The  temperatures  ti,  to,  and  tg  are  the  in¬ 
side  temperature,  outside  air  temperature,  and  the  ground  temperature 
at  basement  depth.  If  the  expression  is  divided  by  (ti-to) ,  a  value 
for  the  transmission  loss  per  hour  per  unit  temperature  difference 
can  be  expressed  as 

Qt/ (ti-to)  =  UgAg  +  UwAw  +  UdAd  +  UcAc  +  { ti-tg) / (ti-to) }UbP  (3) 

As  suggested  in  reference  5,  the  ratio  (ti-tg)/ (ti-to)  is  denoted  by 
Kb  and  assumed  constant,  so  that  (3)  becomes 

Qt/ (ti-to)  =  JgAg  +  UwAw  +  UdAd  +  UcAc  +  Kb  UbP  (4) 

Infiltration  Heat  Loss 

The  heat  load  due  to  infiltrating  air  can  be  determined 
using  the  air  change  method  outlined  in  the  ASHRAE  Handbook  as 

Qi  =  p*Cpa*Vol*I* (Ti-To)  (5) 

p  =  density  of  infiltrating  air 
Cpa  =  specific  heat  of  air,  kJ/ (kg-C) 

Vol  =  building  volume,  meters  cubed 

I  3  number  of  building  air  changes  per  hour 
ti,  to  =  inside  and  outside  air  temperatures,  C 
The  air  change  rate  is  a  function  of  the  density  difference  between 
the  inside  and  outside  air  due  to  temperature  effects  (stack  effect) 
and  the  static  pressure  effect  of  wind  velocity  on  the  structure.  At 
present  there  is  no  analytical  procedure  for  calculating  the  exact 
air  infiltration  into  a  building,  though  numerous  studies  have  been 
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directed  at  this  problem.  Peterson  [23]  has  reviewed  the  results  of 
20  years  of  research  and  produced  an  expression  for  the  number  of  air 
changes  per  hour  in  residential  structures  from  wind  and  stack  effects 
for  different  levels  of  weatherstripping  and  building  construction. 

The  expression  is 

I  =  A  +  B*(ti-to)  +  C*V  (6) 

where  I  is  the  number  of  air  changes  per  hour,  (ti-to)  the  temperature 

difference,  V  the  wind  velocity,  and  A,  B,  and  C  are  empirically  deter¬ 
mined  constants  for  the  level  of  weatherstripping  and  building  tight¬ 
ness.  Peterson’s  expression  was  used  here  for  determining  the  air 
change  rate  and  values  of  the  constants  for  each  of  the  three  house 
types  are  given  in  Table  A-4  in  Appendix  A.  A  value  for  the  number  of 

air  changes  per  hour  for  Model  I  is  calculated  in  the  computer  program 

for  the  building  type  specified  from  the  design  wind  velocity  from 
the  ASHRAE  Handbook  and  the  average  daily  winter  temperature,  both  of 
which  are  input  parameters.  The  air  change  rate  varies  from  about 
0.4  for  a  tightly  constructed  house  (Type  III)  to  1  for  a  loosely  con¬ 
structed  house  (Type  I).  With  the  air  change  rate  thus  specified,  the 
infiltration  heat  loss  can  be  expressed  as  a  loss  per  unit  temperature 
difference  as 

Qi/(ti-to)  =  p*Cpa*Vol*I  (7) 

Building  Heat  Loss  Coefficient 
The  heat  losses  through  transmission  and  infiltration  can 
be  combined  to  give  an  overall  building  heat  loss  per  hour  per  unit 
temperature  difference  as 
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HL  =  (Qt  +  Qi) / (ti-to) 
The  hourly  building  heat  loss  is  then 

(  HL*(ti-to)  ,  for  to<ti 


QLOAD  =  l  (8) 

0  ,  otherwise 

. 

where  QLOAD  is  the  total  heat  loss  in  kilojoules.  Following  the  usual 
practice  [26],  a  reference  temperature  (18  C)  less  than  the  inside  air 
temperature  is  used  to  account  for  the  existence  of  heat  sources  not 
included  in  the  model,  such  as  solar  heat  gain  through  the  windows. 

This  simplified  model  provides  an  approximation  to  the  hourly 
heat  loss  for  the  building.  The  difficulty  with  this  representation 
is  that  it  incorrectly  predicts  the  time  distribution  of  the  heating 
load.  Gutierrez,  et  al.  [9]  found  that  changes  in  the  time  distribu¬ 
tion  of  load  produced  variations  of  as  much  as  17%  in  required  auxil¬ 
iary  energy  in  a  study  of  the  use  of  solar  energy  for  hot  water.  The 
study  considered  only  hot  water  requirements,  and  only  for  a  one  month 
period,  and  so  did  not  address  the  long  term  effects  or  space  heating. 
The  degree  of  departure  of  the  Model  I  approximation  fr^m  the  heat 
loss  given  by  a  more  complex  model  of  space  heating  requirements  will 
depend  primarily  upon  the  magnitude  of  the  hourly  differences  in  infil¬ 
tration  load  from  the  constant  rate  assumed  above,  the  effects  of  so¬ 
lar  radiation  on  the  building,  and  the  variation  of  the  surface  film 
resistance  with  wind.  These  effects  are  considered  in  the  second 
model  of  building  heat  loss  used  in  this  study. 


Model  II 

In  expanding  the  number  of  parameters  in  the  building  heat 


loss  model  the  intent  was  to  include,  in  addition  to  temperature,  two 
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other  major  effects,  wind  and  solar  radiation,  which  act  on  the  build 
ing.  Steady  state  conditions  for  each  hour  of  the  simulation  are 
still  assumed,  so  that  heat  capacity  effects  are  not  included.  The 
degree  to  which  the  inclusion  of  heat  capacity  effects  would  change 
the  results  is  uncertain,  though  McQuiston  and  Parker  [18]  suggest 
they  are  significant  only  in  structures  of  large  thermal  mass,  which 
is  not  the  case  in  this  study. 


Radiation  Effects 

The  effects  of  radiation  on  the  building  were  simulated  by 
the  replacement  of  outside  air  temperature  with  the  sol-air  tempera¬ 
ture  in  computations  of  transmission  heat  loss.  The  sol-air  tempera¬ 
ture,  first  introduced  by  Mackey  and  Wright  [17],  is  a  fictitious 
temperature  that  in  the  absence  of  all  radiative  exchanges  would  pro¬ 
duce  the  same  rate  of  heat  transfer  at  the  exterior  surface  as  actual 
ly  occurs  with  solar  radiation  acting  on  the  surface.  It  can  be  ex¬ 
pressed  as 

Te  =  To  +  al/ho  -  cAR/ho 


where 


Te 

To 

I 

n 

e 

ho 

AR 


sol-air  temperature,  K 
outside  air  temperature,  K 

total  radiation  incident  upon  the  surface,  kJ/(sq  m-h) 

absorptivity  of  the  surface,  dimensionless 

emissivity  of  the  surface,  dimensionless 

surface  film  conductance,  kJ/(sq  m-h-K) 

difference  between  long  wave  radiation  incident  on 
the  surface  and  radiation  from  a  black  body  at 
To,  kJ/ (sq  m-h) . 
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The  last  term  has  been  shown  by  McQuiston  and  Parker  [18]  to  be  neg¬ 
ligible  in  comparison  to  the  first  two  and  in  practice  the  expression 
is  simplified  to 

Te  =  To  +  al/ho. 

The  sol-air  temperature  is  thus  dependent  upon  the  absorptivity  of 
the  surface  and  its  orientation  relative  to  the  sun,  as  well  as  the 
outside  air  temperature.  A  value  of  0.7  was  taken  as  representative 
[24]  of  the  absorptivity  for  the  exterior  building  materials  used  in 
this  study.  The  values  for  outside  surface  film  conductance  were 
those  used  in  the  computation  of  the  overall  heat  transmission  co¬ 
efficient  (U-value)  for  the  surface  under  consideration.  The  solar 
radiation,  I,  incident  upon  the  surface  can  be  calculated  [14]  from 

I  =  RDN*cos0  +  0 . 5*RDU* (1+cosS)  +  0.5*RHT*p*(l-cosB)  (8), 

RDN  =  direct  normal  radiation,  kJ/sq  m 
RDU  =  diffuse  radiation,  kJ/sq  m 
RHT  =  global  horizontal  radiation,  kJ/sq  m 

8  =  tilt  of  the  surface  from  the  horizontal,  deg 
0  =  surface  reflectivity,  dimensionless,  deg 
6  =  angle  of  incidence  of  direct  radiation,  deg. 

The  angle  of  incidence  on  the  surface  is  given  by  [3] 

cos0  =  sin6sin<[>cos8  -  sin6cos<{>sin8cosY  (9) 

+  cos5cost}>cosScosw  +  cosSsin^sinBcosycosw 
+  cos6sin8sinysinti) 

6  =  23 .45*sin[360*(284  +  n)/36o],  deg 

n  =  day  of  the  year 

4>  =  latitude  of  the  station,  deg 
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3  =  tilt  angle  from  the  horizontal,  deg 

y  =  surface  azimuth  angle,  south  =  0,  east  +,  west  -  ,  deg 
u)  =  -  12)/24,  the  hour  angle,  deg 

t  =  solar  time,  h 

Using  (9),  (8),  and  (7)  the  angle  of  incidence,  incident  radiation, 
and  sol-air  temperature  are  calculated  in  the  simulation  each  hour 
the  sun  is  above  the  horizon.  The  transmission  losses  are  then  taken 
as  a  function  of  these  sol-air  temperatures  instead  of  the  outside  air 
temperature. 

Solar  radiation  also  influences  the  hourly  building  heat 
balance  by  the  solar  heat  gain  through  the  windows.  If  a  simplified 
model  of  solar  heat  gain,  with  the  absorptivity  of  the  glass  assumed 
negligible,  is  used,  the  rate  of  heat  entry  from  solar  radiation 
through  the  window  glass  may  be  expressed  [21]  as 

SHG  *  RDN*TAUW*cos(9)  +  0.5*RDU*TAUD  (10) 

where  RDN,  RDU,  and  9  are  as  described  previously,  and  TAUW  and  TAUD 
are  the  transmissivities  of  the  glass  for  direct  and  diffuse  radiation. 
As  suggested  by  Parmelee,  et  al.  [22]  for  vertical  windows,  the  trans¬ 
missivity  for  diffuse  radiation  is  taken  as  a  constant  value  of  0.70 
for  single  pane  and  0.62  for  double  pane  windows.  The  transmissivity 
for  direct  radiation  was  approximated  here  by  a  constant  value  for 
incidence  angles  less  than  60  degrees  and  as  a  linearly  decreasing 
function  of  incidence  angle  for  angles  greater  than  60  degrees.  Fig¬ 
ure  2  depicts  the  transmissivity  of  one,  two,  and  three  glazings  for 
direct  radiation  as  given  by  Duffie  and  Beckman  [6]  (solid  lines)  and 
the  approximations  used  in  this  study  (dashed  lines) .  For  this  model, 


TRANSMISSIVITY 
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Figure  2.  Transmissivity  of  1,2,  and  3  glazings  as 
given  by  Duffie  and  Beckman  (solid  lines)  and  approx¬ 
imations  used  (dashed  lines)  . 
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the  windows  are  assumed  to  have  no  setback  from  the  building  surface, 
and  any  effect  from  drapes  or  curtains  on  either  the  solar  gain  or 
transmission  loss  through  the  windows  is  neglected.  In  addition,  the 
model  assumes  that  window  awnings  which  completely  shade  the  windows 
are  installed  on  June  1  and  remain  in  place  until  September  1,  so  that 
only  the  diffuse  component  is  considered  during  this  period. 

Wind  Effects 

The  effect  of  variations  in  wind  velocity  on  building  heat 
loss  is  considered  in  Model  II  by  including  its  influence  on  surface 
film  conductance  and  on  the  infiltration  rate.  The  convective  com¬ 
ponent  of  the  outside  surface  film  conductance  is  calculated  for  each 
hour  for  each  exposed  surface  as  a  function  of  wind  speed  and  direc¬ 
tion  using  an  expression  developed  by  Kimura  and  reported  in  reference 
20.  In  this  method,  the  air  velocity  close  to  the  surface  is 
VC  =  0.25*V  for  V  >  2  mps 

VC  =  0.5  for  V  <  2  mps 

for  a  relative  wind  direction  less  than  90  degrees,  and 

VC  =  0.3  +  0. 05*V 

for  a  relative  wind  direction  greater  than  or  equal  to  90  degrees, 
where 

V  is  the  reported  wind  velocity,  mps. 

The  relative  wind  is  calculated  from  the  azimuth  angle  of  the  build¬ 
ing  surface  and  the  reported  wind  direction  as 

RWD  =*  ABS  (WAZ  -  DIR)  for  RWD  <  180  degrees 
RWD  -  360  -  RWD  for  RWD  >  180  degrees, 


22 


where  ABS  indicates  the  absolute  value,  RWD  is  the  relative  wind  di¬ 
rection,  WAZ  is  the  surface  azimuth  angle  of  the  building  surface 
(facing  north  =  0),  and  DIR  is  the  reported  wind  direction  in  degrees. 
The  convective  component  of  the  outside  surface  film  conductance,  he, 
is  then  calculated  from 

he  =  67 .049(VC°*605)  (12) 

For  each  hour  of  the  simulation,  the  surface  film  conductance  is  cal¬ 
culated  from  a  convective  component  computed  as  outlined  above  and  a 
constant  value  of  14  kJ/(sq  m-h-K)  for  the  radiation  component  [2]. 

The  heat  transmission  coefficient  (U-value)  for  each  building  surface 
exposed  to  the  wind  is  calculated  from  this  value  of  surface  film  con¬ 
ductance  and  the  same  values  of  thermal  resistance  for  the  interior 
construction  used  in  Model  I.  The  transmission  heat  losses  for  Model 
II  for  each  exposed  surface  are  then  calculated  as  a  function  of  these 
U-values  and  the  sol-air  temperature.  Finally,  infiltration  losses 
for  Model  II  differed  from  Model  I  in  that  an  air  change  rate  was 
calculated  hourly  as  a  function  of  wind  and  temperature  from  equation 
6  and  the  infiltration  heat  loss  calculated  from  equation  7. 

Basement  and  Ceiling  Transmission  Losses 
The  transmission  losses  for  the  ceiling  and  basement  for 
Model  II  also  differ  from  Model  I.  The  transmission  loss  through  the 
ceiling  was  calculated  as  a  function  of  attic  temperature,  calculated 
from  the  following  expression  given  in  the  ASHRAE  Handbook  [2]; 

Tattic  *  (Ac*Uc*Tc+Ta*Ar*Tr+Ta*Aw*Uw)/(Ac*Uc+Ar*Ur+Aw*Uw)  (13) 

where  A  and  U  refer  to  the  areas  and  U-values  and  the  subscripts  c,  r. 
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w  refer  to  the  ceiling,  roof,  and  attic  walls  respectively.  Tattic, 
Tc,  and  Ta  are  the  attic,  ceiling  and  outside  air  temperatures.  As 
suggested  in  reference  2,  the  ceiling  temperature  was  taken  as  1  de¬ 
gree  C  higher  than  the  reference  room  temperature  (20  degrees  C  for 
Model  II) .  The  outside  air  temperature  was  replaced  by  the  applicable 
roof  or  siding  sol-air  value.  Basement  transmission  losses  were  cal¬ 
culated  for  Model  II  by  the  method  outlined  for  Model  I,  except  that 
a  value  of  the  ground  temperature  was  calculated  for  each  day  by 
assuming  it  was  a  sinusoidal  function  of  the  day  of  the  year.  Values 
of  283  K  for  winter  and  294  K  for  summer  from  reference  20  for  average 
values  at  Lake  Hefner,  Oklahoma  were  used  to  determine  the  amplitude. 


Domestic  Hot  Water  Energy  Requirements 
Energy  requirements  for  hot  water  for  both  models  were  cal¬ 
culated  by  assuming  a  requirement  of  75  kg  of  hot  water  per  occupant 
per  day,  evenly  distributed  between  the  hours  of  6  a.m.  and  9  p.m. 

The  number  of  occupants  was  assumed  to  be  4,  5,  and  6  for  house  sizes 
I,  II,  and  III  respectively.  The  required  hot  water  temperature  was 
set  at  60  degrees  C  (140  degrees  F)  and  the  temperature  of  the  feed 
water  from  the  main  assumed  to  be  equal  to  the  ground  temperature. 

The  energy  required  for  hot  water  is  then 

Mw*Cpw(Thw  -  Tg)  ,  6<t<21 


QSHW  = 


(14) 


j  0  ,  otherwise 

QSHW  *  hot  water  energy  requirement,  kJ/h 
Mw  *  mass  of  water,  kg/h 

Cpw  3  specific  heat  of  water,  4.189  kJ/(kg-K) 
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THW  =  hot  water  temperature,  K 
Tg  =  feed  water  temperature,  K 

Summary 

This  chapter  has  presented  the  models  used  for  determining 
the  building  energy  requirements  for  space  and  hot  water  heating. 

With  Model  I  the  hourly  building  heat  requirements  will  vary  only  with 
temperature.  With  Model  II  the  hourly  heat  requirements  will  vary 
with  radiation  and  wind  influences  as  well  as  temperature.  The  solar 
collector  and  storage  models  used  in  the  hourly  simulation  of  system 
performance  will  be  discussed  in  Chapter  V. 


CHAPTER  V 


SOLAR  COLLECTOR  AND  STORAGE  MODELS 

Collector  Model 

The  solar  collector  modeled  is  a  flat  plate  solar  collector 
with  basic  components  as  illustrated  in  Figure  3.  An  absorbing  plate 
is  covered  by  one  or  more  transparent  cover  plates  to  reduce  convec¬ 
tive  heat  losses  and  insulated  to  reduce  transmission  losses  through 
the  back.  The  absorbing  plate,  coated  with  a  material  having  high 
absorptivity  in  the  solar  spectrum,  heats  up  and  in  turn  heats  a  fluid 
flowing  through  the  tubes  in  the  plate.  The  amount  of  energy  collected 
is  a  function  of  the  incident  radiation,  collector  construction,  the 
fluid  temperature,  and  the  surrounding  environmental  conditions. 

In  this  study,  the  collector  is  a  water  heating  flat  plate 
collector  with  the  following  characteristics,  taken  from  Duffie  and 
Beckman's  book  on  solar  energy  thermal  processes  [6]: 

1.  The  absorbing  plate  is  a  copper  tube-in-sheet  plate  with 
10  cm  spacing  between  tubes  and  a  non-selective  coating  which  has  a 
constant  absorptivity,  independent  of  angle  of  incidence,  of  0.95. 

2.  The  glass  covers  (1,  2,  and  3  covers  are  considered)  are 
water  white  glass  with  a  thickness  of  0.3125  cm  and  an  extinction  co¬ 
efficient  of  0.04  per  cm. 
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Figure  3.  Solar  collector  cross  section. 
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3.  The  absorbing  plate  has  a  plate-to-fluid  heat  transfer 
coefficient  of  300  watts  per  sq  meter  per  K. 

Collector  performance  was  modeled  by  the  generalized  flat 
plate  performance  equation  proposed  by  Hottel  and  Whillier  [10]. 

Duffie  and  Beckman  [6]  present  a  complete  derivation  of  this  expres¬ 
sion,  which  relates  the  useful  energy  collected  to  the  collector  and 
environmental  parameters  as  follows 

QU  =  AC*FR* (S-UL*(TCI-TA) )  (15) 

AC  =  collector  area,  sq  m 
FR  *  heat  removal  factor,  dimensionless 
S  =  radiation  absorbed  by  the  plate,  kJ/(sq  m-h) 

UL  =  collector  loss  coefficient,  kJ/(sq  m-h) 

TCI  *  inlet  temperature  of  fluid,  K 
TA  =  ambient  air  temperature,  K. 

The  useful  energy,  QU,  is  the  energy  collected  by  the  heat  transfer 
fluid  In  passing  through  the  tubes  in  the  collector  plate.  The  ra¬ 
diation  absorbed  by  the  plate  is  a  function  of  the  incident  radiation, 
the  transmissivity  of  the  glass  covers,  and  the  absorptivity  of  the 
plate.  The  incident  radiation  and  the  transmissivity  were  calculated 
as  described  in  Chapter  IV.  The  quantities  FR  and  UL  depend  upon  the 
construction  characteristics  of  the  collector,  the  heat  transfer 
characteristics  of  the  plate,  the  mass  flow  rate  through  the  collector, 
and  the  number  of  collector  covers.  Values  of  UL,  the  collector  heat 
loss  coefficient,  for  the  collector  characteristic  previously  listed, 
were  taken  from  graphs  in  Duffie  and  Beckman’s  book  [6]  as 
UL  *  7.2  kJ/(sq  m-K)  for  3  covers 
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UL  =  14.4  kj/(sq  m-K)  for  2  covers 
UL  =  28.8  kj/(sq  m-K)  for  1  cover 

The  effect  of  the  heat  removal  factor,  FR,  in  equation  15  is  to  reduce 
the  calculated  energy  gain  from  the  value  given  by  assuming  the  fluid 
temperature  throughout  the  collector  to  be  at  the  inlet  temperature  to 
its  actual  value,  given  by  a  fluid  temperature  which  increases  as  the 
fluid  flows  through  the  plate.  It  can  be  expressed  [6]  as 

FR  =  (FLO*CP/UL) * (1  -  exp(-UL*F7(FL0*CP))  (16) 

FLO  =  mass  flow  rate  per  unit  area,  kg/sq-m 
CP  =  specific  heat  of  fluid,  kJ/(kg-K) 

UL  =  collector  loss  coefficient,  kJ/(sq  m-K) 

F'  =  collector  efficiency  factor,  dimensionless 
As  demonstrated  by  Duffie  and  Beckman  [6],  F'  is  primarily  a  function 
of  the  physical  characteristics  of  the  absorber  plate  and  can  be  con¬ 
sidered  a  constant.  Values  were  taken  from  Duffie  and  Beckman's  book 
for  the  collector  characteristics  con„:’dered  here  as 
F'  =  0.96  for  1  cover 
F^  =  0.98  for  2  covers 
F^  =  0.99  for  3  covers 

Winn  [24]  reports  that  an  optimum  value  for  the  mass  flow  rate  per 
unit  area  of  collector  can  be  determined  from 

FLO  =■  F'*F'*UL/(2*CP*(F'-FR))  (17) 

Equations  16  and  17,  with  UL,  CP,  and  F"  are  specified,  form  a  set  of 
equations  in  the  two  unknowns,  FLO  and  FR,  which  is  solved  in  the  pro¬ 
gram  using  an  iterative  procedure  for  the  optimal  mass  flow  rate,  FLO, 
and  the  heat  removal  factor,  FR.  With  these  parameters  specified. 
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equation  15  can  be  used  to  determine  the  useful  energy  collected  for 
each  hour  the  solar  system  is  operated. 

The  collectors  modeled  in  this  simulation  were  "operated" 
when  the  temperature  of  the  absorbing  plate  was  heated  to  a  tempera¬ 
ture  ten  degrees  higher  than  the  storage  tank  temperature  (discussed 
below).  As  derived  in  reference  6,  the  temperature  of  the  plate  for 
any  given  hour  can  be  expressed  as 

TP  =  TA  +  S/m  = [S/UL- (TPP-TA) ] *exp (-UL/MCE) 
where  TP  is  the  plate  temperature  for  the  current  hour,  TPP  the  plate 
temperature  for  the  previous  hour,  S,  UL,  TA  as  previously  indicated, 
and  MCE  is  the  effective  heat  capacity  of  the  collector,  calculated 
as  outlined  by  Duffie  and  Beckman  [6]  for  the  1,  2,  and  3  cover  col¬ 
lectors  to  be  6.6,  10.7,  and  15.2  kj/(sq  m-K) . 

Storage  Model 

In  this  study  water  with  a  specific  heat  of  4.189  kJ/(kg-K) 
is  used  as  both  the  heat  transfer  fluid  and  the  storage  medium.  As 
shown  in  L'o'f  and  Tybout's  [16]  study  on  solar  space  heating,  the  op¬ 
timum  mass  of  storage  is  in  the  range  of  50  to  75  kg  per  square  meter 
of  collector,  with  variations  in  this  range  having  a  relatively  small 
effect  on  the  cost  of  delivered  energy  or  the  fraction  of  heat  load 
carried  by  the  solar  system.  The  mass  of  storage  here  is  fixed  at 
61  kg  of  water  per  square  meter  of  collector  as  suggested  by  Butz,  et 
al.  [4].  Beckman,  et  al.  [3]  show  that  modeling  stratification  effects 
in  water  storage  produces  results  only  1  to  3%  different  from  model¬ 
ing  an  unstratified  system,  so  the  system  modeled  here  uses  an  un¬ 
stratified  model  for  storage. 


30 


The  storage  vessel  is  a  cylindrical  water  tank  taken  to  have 
a  height  to  diameter  ratio  of  3:1  and  insulation  with  a  heat  trans¬ 
mission  coefficient  of  1.44  kJ/(sq  m-hr-K)  [6].  The  size  and  surface 
area  are  determined  in  the  program  from  the  mass  of  storage  water  and 
the  assumptions  of  a  constant  density  of  1000  kg  per  cubic  meter  and 
the  height  to  diameter  ratio  of  3:1.  Thermal  losses  from  storage  are 
then  determined  from  the  area  U-value  product,  the  storage  temperature 
and  the  room  temperature. 

A  schematic  of  the  solar  collector  and  storage  system  illus¬ 
trating  the  mode  of  operation  is  presented  in  Figure  4.  Energy  is 
transferred  from  storage  to  incoming  water  from  the  main  for  hot  water 
use  by  means  of  a  heating  coil  located  within  the  solar  system  storage 
tank,  and  from  storage  to  the  house  by  means  of  a  water  to  air  heat 
exchanger.  The  size  and  heat  transfer  characteristics  of  the  heat 
exchanger  are  not  specifically  modeled,  but  it  is  assumed  as  suggest¬ 
ed  in  reference  3,  that  the  water  and  air  flow  rates  are  modulated  so 
that  the  average  rate  of  energy  transfer  across  the  load  heat  exchang¬ 
er  is  equal  to  the  average  rate  of  energy  required  by  the  heating  load 
Thus  the  heat  exchanger  effects  may  be  neglected. 

For  the  system  described  above,  an  energy  balance  on  the 
fully  mixed  storage  tank  is  given  by 

M*Cps*dTs/dT  =  QU  -  QL  -  QSHW  -  QTL  (18) 

M,  Cps  =  mass,  specific  heat  of  storage,  kg,  kJ/(kg-K) 

Ts  =  storage  temperature,  K 
T  =  t ime ,  h 

QU  »  rate  of  transfer  of  energy  to  storage,  kJ/h 
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Figure  4.  Solar  system  schematic. 
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QL  =  rate  of  transfer  of  energy  to  space  heat  load,  kJ/h 
QSHW  =  rate  of  transfer  of  energy  to  hot  water  load,  kJ/h 

QTL  =  rate  of  heat  loss  from  storage,  kj/hr,  given  by 
QTL  =  Us*As*(Ts-Tr) 

Us,  As  =  heat  transfer  coefficient  and  area  of  storage  tank 
Ts,  Tr  =  storage  tank  and  room  temperatures,  K 
As  indicated  by  Duffle  and  Beckman  [6],  equation  18  may  be  integrated 
in  one  hour  time  steps  by  Euler  integration  to  determine  the  storage 
tank  temperature  at  the  end  of  the  hour  from  the  temperature  at  the 
beginning  of  the  hour  and  the  known  energy  inputs  and  outputs.  The 
expression  is  solved  hourly  in  the  simulation  for  the  storage  tempera¬ 
ture,  which  has  an  upper  limit  of  90  degrees  C.  When  the  useful  energy 
collected  raises  the  storage  temperature  above  this  maximum,  the  stor¬ 
age  temperature  is  set  equal  to  the  maximum  value  and  the  amount  of 
excess  energy,  which  is  assumed  dumped  to  the  atmosphere,  tabulated. 

The  simulation  thus  assumes  a  controller  which  senses  the  storage  tem¬ 
perature  and  halts  the  collector  pump  when  the  maximum  temperature  is 
reached.  As  indicated  in  reference  3,  the  minimum  storage  temperature 
at  which  the  water  to  air  heat  exchanger  is  effective  is  on  the  order 
of  26-30  degrees  C.  A  temperature  of  29  degrees  C  was  used  here  as 
the  storage  temperature  below  which  delivery  of  energy  from  storage 
to  the  loads  is  halted  and  auxiliary  energy  used  for  both  space  and 
hot  water.  Thus,  for  space  heating,  either  all  or  none  of  the  energy 
requirements  are  met  from  storage  in  any  given  hour,  based  upon  the 
calculated  storage  temperature  for  that  hour.  The  solar  system  stor¬ 
age  acts  as  a  preheater  for  water  from  the  mains  when  the  storage 
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temperature  is  below  the  required  hot  water  value  of  60  degrees  C. 

If  the  storage  temperature  is  above  60  degrees,  a  control  system  and 
a  mixing  valve,  not  specifically  modeled,  are  assumed  to  operate  to 
mix  water  from  the  mains  with  water  heated  by  the  storage  system  to 
deliver  60  degree  C  water  to  the  domestic  hot  water  tank.  Thus,  for 
domestic  hot  water,  auxiliary  energy  is  used  exclusively  if  the  stor¬ 
age  temperature  is  below  the  minimum  storage  temperature,  is  used  to 
increase  the  domestic  hot  water  temperature  to  the  required  value  of 
60  degrees  if  the  storage  temperature  is  above  minimum  but  below  60 
degrees,  and  is  not  required  if  the  storage  temperature  is  above  60 
degrees. 

With  the  known  requirements  for  space  heat  and  hot  water, 
the  solar  energy  input,  and  the  simulation  of  the  collector  and  stor¬ 
age  system  as  outlined  above,  the  annual  requirements  for  auxiliary 
energy  can  be  determined  and  used  to  find  the  most  economic  mix  of 
conventional  and  solar  energy.  This  economic  optimization  of  the  so¬ 
lar  system  size  is  discussed  in  the  next  chapter. 


CHAPTER  VI 

COST  OPTIMIZATION  OF  COLLECTOR  AREA 

Introduction 

The  preceding  chapters  have  outlined  the  methods  of  model¬ 
ing  the  thermal  performance  of  the  components  used  in  the  solar  sys¬ 
tem  simulation  and  have  not  dealt  with  questions  of  cost.  However, 
the  decision  to  incorporate  solar  energy  for  the  heating  of  a  resi¬ 
dence  and  the  selection  of  the  size  of  the  system,  assuming  the  in¬ 
dividual  thinks  rationally,  is  an  economic  one.  This  chapter  will 
present  the  methodology  used  to  determine  the  optimum  size  of  the  so¬ 
lar  system  from  the  results  of  the  thermal  performance  simulation  of 
the  house  and  solar  system. 

As  previously  indicated  in  Chapter  I,  studies  by  Lof  and 
Tybout  [15,  16]  have  clearly  shown  that  designing  a  solar  system  to 
provide  100%  of  the  heating  requirements  for  a  conventionally  con¬ 
structed  residence  is  not  feasible.  The  size  of  the  collector  array 
and  storage  system  required  raise  the  cost  far  above  the  cost  of  con¬ 
ventional  energy  sources  and  produce  architecturally  undesirable  con¬ 
sequences  [1],  such  as  the  inability  to  incorporate  the  solar  system 
into  the  construction  of  the  residence.  Consequently,  a  conventional 
system  capable  of  meeting  the  worst  case  heating  requirements  of  the 
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building  must  be  part  of  the  residence  and  its  costs  are  required 
whether  or  not  a  solar  system  is  included.  The  only  savings  then 
offered  by  solar  energy,  where  it  is  of  low  enough  cost,  is  in  the 
reduction  of  the  cost  required  for  fuel  for  the  conventional  heating 
system.  Thus,  as  demonstrated  by  Duffie,  Beckman,  and  Klein  in  ref¬ 
erence  3,  the  conventional  system  cost  can  be  considered  as  a  base 
cost  common  to  both  systems  and  the  cost  analysis  carried  out  as  a 
comparison  of  the  costs  above  base. 

In  a  similar  manner,  the  cost  of  insulation  was  included  in 
this  study  as  a  cost  above  base  by  assuming  that  a  decision  regarding 
the  level  of  insulation  is  made  concurrently  with  the  decision  regard¬ 
ing  the  solar  system.  Type  II  insulation,  representing  a  moderately 
insulated  house,  is  taken  as  the  base  level.  The  cost  of  increasing 
the  insulation  to  Type  III,  or  the  savings  from  deciding  to  insulate 
only  to  the  level  of  Type  I  are  calculated  and  included  in  the  cost 
analysis.  In  addition,  it  is  assumed  here  that  the  solar  system  is 
well  designed  and  has  negligible  maintenance  and  pumping  costs. 

Annual  Energy  Cost 

The  annual  cost  of  energy  with  the  mix  of  solar  and  conven¬ 
tional  systems,  under  the  assumptions  listed  above,  may  be  expressed 
as : 

EUAC  =  (SOLCS  +  CSTI) *1  +  CSTA  (19) 

EUAC  s  uniform  annual  cost  of  energy  for  space  and  hot  water, 
dollars  per  year 

SOLCS  *  total  capital  cost  of  solar  system,  dollars 

CSTI  =  total  capital  cost  of  insulation,  above  or  below 
base,  dollars 
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I  =  capital  recovery  factor  for  assumed  interest  rate 
and  time  period,  fraction  per  year 

CSTA  =  uniform  annual  cost  of  auxiliary  energy,  where  auxil¬ 
iary  energy  refers  to  the  energy  provided  by  the  con¬ 
ventional  heating  system. 

Solar  System  Cost 

The  total  capital  cost  of  the  solar  system  is  the  sum  of  the 
costs  for  the  collector,  piping,  pumps,  controls,  and  storage  minus 
the  tax  incentive  provided  under  current  law  for  installation  of  a 
solar  system.  A  cost  of  $100  per  square  meter  for  a  two  cover  collec¬ 
tor  was  taken  as  a  base  collector  cost  [18]  ,  with  an  additional  cost 
or  savings  of  $8.60  for  a  three  cover  and  one  cover  collector,  respec¬ 
tively.  As  indicated  in  the  study  by  Butz,  et  al.  [4],  the  cost  of 
storage  can  be  considered  a  linear  function  of  collector  area  and 
their  figure  of  $8/sq  m  for  the  cost  of  the  storage  vessel,  water, 
and  insulation  was  used  here.  The  additional  solar  system  costs  for 
piping,  pumps,  and  controls,  as  suggested  in  reference  3,  were  set  at 
$375,  independent  of  collector  area.  Current  tax  policy  [25]  provides 
that  of  the  total  cost  of  a  solar  installation,  30%  of  the  first  $2000 
and  20%  of  the  remaining  cost,  up  to  a  maximum  of  $10000,  are  rebated 
to  the  purchaser  as  tax  credit.  Under  these  assumptions,  the  initial 
capital  cost  of  the  solar  system  is 

CCOST  =  ($100  +  $8. 6*(N-2)  +  $8)*AC  +  $375 
where  N  is  the  number  of  collector  covers  and  AC  is  the  collector 
area.  The  tax  savings  is  then 
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TAXSV  = 


0. 3*CC0ST 

$600  +  0. 2* (CCOST-$2000) , 
$2200 


CCOST<$2000 
$2000<CCOST<$ 10000 
CCOST>$10000 


The  final  capital  cost  of  the  solar  system  is  then 


SOLCS  =  CCOST  -  TAXSV 


Insulation  Cost 

As  indicated  previously,  the  cost  of  insulation,  CSTI,  was 
treated  as  a  cost  above  or  below  the  base  cost  of  insulating  the  resi¬ 
dence  to  insulation  Type  II.  The  cost  of  fibrous  insulation  used  for 
the  wall  and  ceiling  insulation  in  the  residences  considered  here  was 
taken  from  the  February  1980  ASHRAE  Journal  [13]  as  $0.08  per  sq  ft 
(0.0929  sq  m)  for  the  first  inch  (2.54  cm)  and  $0.05  per  sq  ft  for 
each  additional  inch  of  insulation.  These  costs,  the  total  wall  and 
ceiling  area,  and  the  difference  in  insulation  thickness  from  Type  II 
insulation  were  then  used  to  calculate  the  total  cost  or  savings  above 
or  below  base  for  each  insulation  type  for  each  house  size  considered. 
Table  1  summarizes  the  results  obtained.  The  capital  recovery  factor, 
I,  used  to  convert  the  total  capital  costs  to  a  uniform  annual  cost 
was  0.11746  [8]  for  the  assumed  20  year  amortization  period  at  an 
interest  rate  of  10%. 


Auxiliary  Energy  Cost 

The  cost  of  auxiliary  energy  is  determined  from  the  amount 
of  auxiliary  energy  required  and  the  unit  cost  of  the  fuel  used.  An 
initial  cost  of  $6  per  million  kJ,  at  the  upper  range  of  values  con¬ 
sidered  by  Butz,  et  al.  in  their  1974  study,  was  used.  It  was  not. 


INSULATION  TYPE 
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i 


HOUSE  SIZE 


I 

II 

III 

I 

-  $263 

-  $373 

-  $491 

II 

$0 

$0 

$0 

III 

$717 

$1021 

$1346 

Table  1  .  Capital  cost  of  insulation  above  or 
below  insulation  Type  II.  Negative  values  are 
savings,  positive  values  increased  cost. 
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however,  held  constant  over  the  20  year  period  considered  in  each  sim¬ 
ulation.  A  Westinghouse  study  prepared  for  the  Department  of  Energy 
[28]  indicates  that  conventional  energy  costs  will  increase  at  about 
8%  per  year  through  the  year  2000,  and  this  factor  was  taken  into 
account  in  order  to  make  an  equitable  comparison  between  the  uniform 
annual  cost  of  the  solar  system  and  the  annual  cost  of  auxiliary 
energy  over  the  20  year  period.  This  was  accomplished  by  assuming  the 
value  of  required  auxiliary  energy  obtained  from  the  simulation  using 
the  average  year  represents  the  average  annual  value  over  the  20  year 
period.  Using  the  figure  of  $6  per  million  kJ  as  an  initial  cost, 
the  total  cost  of  auxiliary  energy  for  the  20  year  period  was  calculat¬ 
ed  and  converted  to  an  equivalent  uniform  annual  cost  by  the  use  of  a 
capital  recovery  factor  (0,10185)  at  8%  for  20  years. 

Summary 

Under  the  assumptions  listed  above,  equation  19  can  be  seen 
to  represent  the  total  annual  cost  above  base  for  delivered  energy 
for  space  and  hot  water  heating  from  the  solar  and  conventional  systems 
as  a  function  of  collector  area.  Figure  5  illustrates  the  relative 
contributions  of  the  two  terms  in  the  equation  with  increasing  collec¬ 
tor  area.  As  the  collector  area  increases,  collector  costs  (for  a 
fixed  insulation  cost)  increase  linearly.  However,  as  annual  solar 
system  efficiency  decreases  with  increasing  collector  size,  less  and 
less  auxiliary  energy  is  replaced  by  a  unit  increase  in  collector  area 
because  the  solar  system  does  not  operate  at  full  capacity  on  a  year 
round  basis.  The  result  is  a  minimum  on  the  curve  depicting  the  total 
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Figure  5.  Relative  contributions  of  collector  and 
insulation  and  auxiliary  energy  to  cost  of  delivered 
energy  for  space  and  hot  water  heating. 
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annual  cost  of  energy  which  occurs  at  the  optimum  collector  size. 
Optimum  in  this  case  refers  to  the  collector  size  which  produces  the 
lowest  annual  energy  cost  for  space  and  hot  water  heating  from  the  mix 
of  solar  and  conventional  energy  sources. 

This  method  of  economic  analysis,  coupled  with  the  thermal 
analyses  of  solar  system  performance  and  building  heat  loss  as  out¬ 
lined  in  previous  chapters,  forms  the  basis  for  comparison  of  the 
effect  of  the  two  methods  of  modeling  heat  loss  on  the  selection  of 
optimum  collector  size.  The  results  are  presented  in  the  next  chapter. 


CHAPTER  VII 

RESULTS,  CONCLUSIONS,  AND  RECOMMENDATIONS 

Results 

The  procedures  outlined  in  previous  chapters  were  used  to 
develop  a  numerical  model  for  the  simulation  of  building  heat  loss 
and  solar  system  performance.  The  program,  listed  in  Appendix  C, 
uses  the  Oklahoma  City  average  year  discussed  in  Chapter  I  for  hourly 
simulation  of  the  building  heat  loss  and  solar  system  performance 
using  Model  I  and  Model  II  methods.  It  determines  the  amount  of  aux¬ 
iliary  energy  required,  the  percentage  of  the  heating  requirements 
provided  by  the  solar  system,  and  the  annual  cost  of  energy  for  space 
and  hot  water  heating.  A  total  of  172  simulation  runs  for  various 
combinations  of  house  size,  insulation  type,  and  number  of  collector 
covers  and  collector  area  were  accomplished  and  are  summarized  in 
Tables  B-l  to  B-27  in  Appendix  B. 

Each  table  in  Appendix  B  presents  a  comparison  of  the  simu¬ 
lation  results  for  collector  areas  up  to  100  square  meters  for  a  given 
house  size,  insulation  type,  and  number  of  collector  covers.  Table  2 
is  an  example  of  the  information  contained  in  each  table  in  Appendix 
B.  The  first  column,  labeled  AC,  gives  the  collector  size  in  square 
meters  for  which  the  comparison  between  Model  I  and  Model  II  results 
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AC 

<sq 

m) 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 
ENERGY  j 

(million  | 
k  1 1  o  i  cu  1  e  s) 

15 

I 

$963 

EO 

58.4 

.414 

$237 

$1995 
-  599 

65.4  . 

II 

$947 

23.5 

61.1 

.378 

$253 

$1396 

64.1 

30 

I 

$878 

45.6 

70.3 

.354 

$318 

$3615 
-  923 
$2692 

46.0 

II 

$859 

sa 

72.5 

.326 

$342 

44.4 

45 

I 

$850 

64.2 

76.8 

.328 

$34  3 

$5235 

-1247 

$3988 

31.3 

II 

$835 

MM 

79.5 

.284 

$365 

30.0 

60 

I 

$870 

77.8 

81.8 

.  308 

$318 

$6855 

-1571 

$5284 

20.4 

II 

$865 

n 

83.8 

.261 

$333 

20.0 

75 

I 

$929 

87.2 

85.6 

.291 

$256 

$8475 

-1895 

$6580 

12.7 

II 

$925 

87.4 

87.3 

.244 

$270 

12.5 

ICO 

I 

$1118 

199 

91.0 

.268 

$60 

$11175 
-  2200 
$  8975 

5.2 

II 

$1116 

96.4 

92.2 

.223 

$71 

5.1 

Table  2.  Example  of  tables  contained  in  Appendix  B.  This  table 
summarizes  the  simulation  results  for  house  size  III,  insulation 
type  II,  2  cover  collector. 
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is  made.  The  next  column  lists  the  uniform  annual  cost  of  energy  for 
space  and  hot  water  requirements  (EUAC,  equation  19)  obtained  using 
each  model  for  that  collector  size.  In  the  next  two  columns,  the 
percentage  of  the  total  annual  energy  requirement  for  space  heat  and 
hot  water  provided  by  the  solar  system  are  listed  for  each  model. 

Next  is  the  annual  efficiency  of  the  solar  system,  defined  as  the 
total  amount  of  useful  energy  provided  by  the  solar  system  divided  by 
the  annual  amount  of  energy  incident  upon  the  collector  array.  The 
column  labeled  Annual  Savings  gives  the  dollar  value  of  the  annual 
energy  saved  by  the  solar  system  for  the  given  collector  area  above 
that  given  by  the  use  of  conventional  fuel  alone  for  space  and  hot 
water  heating.  Under  the  heading  Solar  System  Cost,  the  figures  rep¬ 
resent  the  initial  capital  cost  of  the  solar  system,  minus  the  tax 
credit  provided,  and  the  resulting  figure  for  the  final  capital  cost 
of  the  solar  system  (SOLCS,  equation  19).  Finally,  the  last  column 
presents  the  annual  amount  of  auxiliary  energy  which  the  simulation 
produces  for  each  heat  loss  model  in  millions  of  kilojoules. 

It  is  apparent,  upon  reviewing  the  results  given  in  these 
tables,  that  the  same  optimum  collector  area  (as  defined  in  Chapter 
VI)  would  be  selected  for  each  building  heat  loss  model  with  a  given 
house  size,  insulation  level,  and  number  of  collector  covers.  The 
results  from  Model  II  indicate  a  slightly  lower  annual  energy  cost 
and  amount  of  auxiliary  energy  and  a  slightly  higher  annual  savings 
than  the  Model  I  results,  though  the  difference  is  generally  less 
than  5Z.  Model  II  results  also  show  a  slightly  greater  percentage  of 
the  space  and  hot  water  heating  loads  carried  by  the  solar  system 
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than  indicated  by  Model  I  results,  but  the  annual  system  efficiency 
for  Model  II  is  lower  than  with  Model  I,  a  point  which  will  be  dis¬ 
cussed  in  more  detail  later.  It  is  also  clear  that  regardless  of  the 
heat  loss  model  used  or  the  size  of  the  house  and  solar  system,  a  two 
cover  collector  is  always  significantly  more  cost  effective  than  one 
cover  and  slightly  more  cost  effective  than  three  covers,  a  result 
which  Lof  and  Tybout  found  to  be  true  in  their  study  of  solar  space 
heating  [15].  Finally,  recalling  that  the  figure  given  for  the  annual 
cost  of  energy  includes  the  cost  of  insulation  above  or  below  base, 
the  results  show  that  both  models  predict  that  the  highest  insulation 
level  provides  the  least  cost  of  energy  for  space  and  hot  water  heat¬ 
ing. 

These  results  are  summarized  in  Figures  6,  7,  and  8,  which 
present  plots  of  the  annual  cost  of  energy  for  space  and  hot  water 
heating  as  a  function  of  collector  area  for  house  sizes  I,  II,  and 
III,  respectively.  The  label  on  the  cost  curves  in  each  plot  indi¬ 
cate  the  insulation  level,  I,  II,  or  III  and  the  number  of  collector 
covers,  1  or  2.  Since  the  values  for  both  heat  loss  models  resulted 
in  nearly  identical  plot  curves,  for  clarity  the  curves  for  Model  II 
only  are  presented.  Similarly,  the  values  for  a  3  cover  collector 

were  nearly  coincident  with  the  two  cover  values  in  each  case  and 

were  not  plotted.  These  figures  clearly  demonstrate  the  cost  effec¬ 
tiveness  of  increased  insulation  for  each  house  size,  and  the  superi¬ 
ority  of  the  two  cover  collector  over  the  one  cover  collector. 

Figures  9,  10,  and  11  are  presented  as  an  example  of  the 

daily  performance  of  the  two  models  of  heat  loss  and  the  influence  bv 


ANNUAL  ENERGY  COST, DOLLARS 


46 


o 


Figure  6.  Annual  cost  of  delivered  energy  for 
space  and  hot  water  heating  for  house  size  I  for 
3  insulation  levels  and  1  or  2  collector  covers. 
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Figure  7.  Annual  cost  of  delivered  energy  for 
space  and  hot  water  heating  for  house  size  II  for 
3  insulation  levels  and  1  or  2  collector  covers. 
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COLLECTOR  AREA  (SQ  M) 

Figure  8.  Annual  cost  of  delivered  energy  for 
space  and  hot  water  heating  for  house  size  III  for 
3  insulation  levels  and  1  or  2  collector  covers. 
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each  on  the  solar  system  performance.  They  were  obtained  with  house 
size  II,  insulation  type  II,  a  two  cover  collector,  and  36  square 
meters  of  collector.  Figure  9  shows  the  daily  average  value  of  tem¬ 
perature  in  degrees  C,  the  daily  average  wind  speed  in  meters  per 
second,  and  the  daily  average  cloud  cover  in  tenths.  In  Figure  10, 
the  upper  two  plots  are  the  daily  average  values  of  building  heat 
requirements  for  space  and  hot  water  and  the  lower  two  plots  are  the 
daily  average  values  of  auxiliary  energy,  labeled  for  I  for  Model  I 
and  II  for  Model  II.  The  values  of  building  heat  required  and  auxil¬ 
iary  energy  are  nearly  coincident  for  both  models,  with  the  greatest 
departures  in  the  spring  and  fall  when  the  building  space  heat  demand 
is  at  a  low  level.  Figure  11  shows  the  daily  average  storage  tempera¬ 
tures  for  both  models.  It  indicates  that  the  storage  temperature  with 
Model  II  is  generally  higher  than  with  Model  I.  The  lower  annual  so¬ 
lar  system  efficiency  referred  to  earlier  is  a  direct  consequence  of 
this  higher  storage  temperature,  since  the  collector  losses  are  great¬ 
er  and  the  useful  energy  collected  lower  with  a  higher  fluid  inlet 
temperature  to  the  collector. 

A  more  detailed  picture  is  given  in  Figures  12,  13,  14,  and 
15  where  a  48  hour  period  of  the  hourly  simulation  for  each  model  is 
shown  for  the  same  house  and  solar  system.  Figure  12  is  an  hourly 
plot  of  the  same  weather  parameters  as  given  in  Figure  9.  In  Figure 
13,  the  hourly  heat  requirements  using  Model  I  and  Model  II  are  pre¬ 
sented,  along  with  the  auxiliary  energy  required  for  both  models  and 
the  solar  gain  through  the  building  windows  for  Model  II.  The  Model 
II  values  of  heat  loss  are  greater  than  Model  I  at  night  and  in  the 


TIME,  HOURS 

Figure  12.  Hourly  values  of  temperature  in  degrees  C,  wind  speed  in 
meters  per  second,  and  cloud  cover  in  tenths. 
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heat  loss  models  I  and  II  and  useful  solar  gain  for  model  II 
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Figure  14.  Hourly  solar  system  storage  temperature  with  heat  loss  models  I  and  II 
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early  morning,  but  are  significantly  reduced  below  the  Model  I  values 
by  the  solar  gain  during  the  day.  The  effect  on  the  storage  tempera¬ 
ture  and  efficiency  of  the  solar  system  is  apparent  in  Figures  14  and 
15.  As  seen  in  Figure  14,  the  storage  temperature  with  Model  II  fall 
below  the  value  given  by  Model  I  at  night  and  rises  above  the  value 
given  by  Model  I  during  the  day  when  solar  energy  collection  and  in¬ 
put  to  storage  takes  place.  The  effect  of  the  solar  system  efficien¬ 
cy  can  be  seen  in  Figure  15  where  the  total  energy  incident  on  the 
collector,  the  useful  energy  collected,  and  the  collector  losses  are 
plotted  for  each  heat  loss  model.  With  the  same  value  of  incident 
energy  on  both  collectors,  the  collector  losses  with  Model  II  exceed 
those  with  Model  I  because  of  the  higher  fluid  inlet  temperature. 
Consequently,  less  useful  energy  is  collected  and  the  system  efficien 
cv  is  lower. 


Conclusions 

These  results  show  that  the  time  distribution  of  heating 
demand  as  simulated  by  heat  loss  Model  II  does  result  in  a  signifi¬ 
cant  difference  in  solar  system  performance  on  a  short  term  basis 
from  the  results  given  by  Model  I.  However,  in  contrast  to  the  re¬ 
sults  found  by  Gutierrez,  et  al.  [9]  for  domestic  hot  water  only, 
these  short  term  differences  between  the  two  models  do  not  result  in 
significantly  different  long  term  fractions  of  the  load  carried  by 
the  solar  system.  Indeed,  it  is  the  departure  in  the  time  distribu¬ 
tion  of  the  heat  load  given  by  Model  II  from  that  of  Model  I  which 


produces  a  lower  annual  efficiency  and  results  in  the  selection  of 
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the  same  optimum  collector  size  for  both  approaches  even  though  the 
Model  II  approach  gives  lower  annual  auxiliary  requirements.  Thus, 
the  Model  II  approach,  though  more  complex  and  requiring  much  more 
computer  time,  does  not  produce  different  results  in  the  selection  of 
the  solar  system  size  at  Oklahoma  City.  This  indicates  that  where 
the  purpose  is  the  simulation  of  long  term  average  results  for  the 
selection  of  the  size  of  the  solar  system,  the  increased  complexity 
and  cost  of  the  Model  II  approach  are  not  worthwhile.  On  the  other 
hand,  if  the  purpose  is  the  simulation  of  the  short  term  performance 
characteristics  of  a  solar  system,  these  results  indicate  that  a  high¬ 
er  level  approach  than  that  of  Model  I  is  necessary. 

Recommendations  for  Further  Research 

1.  This  analysis  was  accomplished  at  a  single  station. 
SOLMET  data  tapes  are  available  for  locations  with  significantly  dif¬ 
ferent  climate  regimes  and  could  be  used  to  determine  whether  differ¬ 
ent  climates  would  produce  different  results. 

2.  The  building  models  used  were  of  conventional  construc¬ 
tion  characteristics.  An  extension  of  the  analysis  to  houses  designed 
to  be  solar  residences  could  be  accomplished.  This  could  include  an 
analysis  of  the  effect  of  varying  the  size  and  location  of  the  build¬ 
ing  windows  in  a  conventionally  constructed  residence. 

3.  Each  house  considered  here  was  a  south  facing  structure 
with  the  collector  optimally  tilted  for  space  heating.  The  program 
is  written  to  accept  off-south  orientations  and  any  collector  tilt 
angle  and  could  be  used  to  study  the  differences  between  the  optimum 


south  facing  system  and  other  orientations. 
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Table  A-l.  House  Design  Data 


I 

House  Size 

II 

III 

Area 

92.9  sq  m 

134.8  sq  m 

180.0  sq  m 

Configuration 

7.62  x  12.2  m 

9.5  x  14.3  m 

12.0  x  15.0  m 

Foundation  Length 

39.6  m 

47.5  m 

54.0  m 

Window  Area 

9.7  sq  m 

13.5  sq  m 

17.7  sq  ra 
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Table  A-2.  House  Construction  Data 


Insulation  Type 

I 

II 

III 

Roof 

tile  shingles 
building  paper 
on  wood  sheathing 

tile  shingles 
building  paper 
on  wood  sheathing 

asphalt  shingles 
on  wood  sheathing 

Attic- 

Walls 

lapped  wood  siding 
wood  sheathing 

brick  veneer 
wood  sheathing 

1 

brick  veneer 
wood  sheathing 
in  insulation  board 

Ceiling 

3/8  in  gypsum 
1/2  in  plaster 

1  in  fibrous  insul 

3/8  in  gypsum 

1/2  in  plaster 

3  in  fibrous  insul 

3/8  in  gypsum 

1/2  in  plaster 

12  in  fibrous  insul 

Walls 

gypsum  lath 

4  in  air  space 
building  paper 
wood  sheathing 
lapped  wood  siding 

gypsum  lath 

3  in  air  space 

1  in  fibrous  insul 
building  paper 
wood  sheathing 

4  in  face  brick 

gypsum  lath 

1  in  air  space 

3  in  fibrous  insul 
building  paper 
wood  sheathing 

8  in  common  brick 

Glass 

single  pane 
aluminum  sash 

double  pane 
wood  sash 

double  pane 
wood  sash 

storm  sash  in  winter 

Doors 

standard  1  in 
exposed 

standard  1  1/2  in 
storm  door 

2  in  insulated 
storm  door 

Basement 

below  grade 
concrete 

uninsulated  walls 

below  grade 
concrete 

1  in  wall  insul 

below  grade 
concrete 

2  in  wall  insul 

65 


Table  A-3.  Heat  Transmission  Coefficients,  kJ/(sq  m-h-C) 

Insulation  Type 


I 

II 

III 

Roof 

12.4 

12.4 

10.0 

Attic  Walls 

8.9 

8.9 

5.9 

Ceiling 

3.9 

1.6 

0.9 

Walls 

5.1 

2.8 

1.5 

Glass 

25.4 

10.7 

9.8 

Doors 

7.6 

6.5 

5.7 

Basement 

7.2 

2.8 

2.2 

Table  A-A.  Coefficients  for  Peterson's  Infiltration  Equation 


Insulation 

Coefficient 

Type 

A 

B 

C 

I 

0.10 

0.216 

0.0300 

II 

0.10 

0.0162 

0.0225 

III 


0.10 


0.0108 


0.0150 


APPENDIX  B 

SUMMARY  OF  SIMULATION  RESULTS 
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AC 

(sq 

m) 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 
ENERGY 
(million 
k  i  1  o  3  ou  1  e  s) 

g 

I 

$975 

9.9 

EB 

.357 

B 

• 

■ 

73.4 

II 

$933 

11.5 

49.2 

.306 

$136 

70.0 

I 

$908 

23.3 

56.3 

.278 

$186 

$2164 
-  633 

61.8 

lo 

II 

$886 

25.1 

58.3 

.267 

$192 

$1531 

60.0 

27 

I 

$890 

33.6 

60.2 

.248 

$202 

$3059 

53.4 

II 

$873 

35.9 

63.9 

.  236 

$211 

oil 

$2248 

52.1 

I 

$883 

43.1 

65.0 

.229 

$208 

$3953 
-  991 
$2962 

46.0 

o  b 

II 

$880 

44.6 

67.8 

.209 

$210 

45.7 

50 

I 

$901 

m 

67.7 

.202 

$188 

$5345 

-1269 

$4076 

36.7 

II 

$898 

56.7 

68.3 

.183 

$196 

36.5 

fl 

I 

$982 

71.6 

75.0 

.161 

$103 

$7830 

-1766 

$6064 

24.3 

II 

$980 

72.3 

77.0 

.151 

$113 

24.1 

100 

I 

$1118 

83.1 

80.0 

.131 

$  0 

$10315 
-  2200 
$  8115 

15.7 

II 

$1115 

83.3 

80.7 

.130 

$  0 

15.1 

Table  B-l.  Summary  of  simulation  results  for  house  size  I,  insulation 
type  I,  1  cover  collector. 
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AC 

vsq 

m) 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

(million 

kilojoules) 

Q 

I 

$937 

13.6 

52.5 

.433 

$158 

$1347 
-  404 
$  942 

69.8 

II 

$893 

n 

54.7 

.396 

$176 

66.2 

18 

! 

I 

$847 

30.6 

63.9 

.375 

$247 

$2319 

aaA 

55.6 

II 

$819 

32.5 

64.5 

.3  54 

$259 

$1655 

53.3 

1 

I 

$810 

43.8 

69.8 

.351 

$281 

$3291 

45.1 

II 

$781 

46.1 

71.4 

.324 

$301 

42.8 

36 

I 

$788 

55.7 

74.3 

.335 

$301 

$4263 

-1053 

$3210 

35.8 

II 

$761 

58.0 

76.7 

.297 

$324 

33.6 

!  50 

I 

$791 

70.5 

79.6 

.315 

$294 

$5775 

-1355 

$4420 

24.4 

II 

$772 

72.1 

81.7 

.267 

$315 

22.9 

75 

I 

$879 

88.0 

86.2 

.287 

$198 

$8475 

-1895 

$6580 

10.9 

II 

$871 

88.7 

87.8 

.  240 

$215 

10.2 

i  nr\ 

1 

$1074 

96.8 

91.6 

.264 

$  o 

3.8 

1UU 

11 

m 

96.7 

92.8 

.220 

$  0 

3.7 

Table  b-2.  Summary  of  simulation  results  for  house  size  I,  insulation 
type  I,  2  cover  collector. 


AC 

m) 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

(million 

kilojoule^ 

9 

I 

$943 

13.7 

52.6 

.459 

$152 

$1424 
-  427 
$  997 

69.7 

II 

$895 

15.7 

55.6 

.414 

$173 

65.9 

18 

I 

$857 

30.9 

64.2 

.430 

$235 

$2474 
-  695 
$1779 

55.5 

II 

$820 

33.2 

67.0 

.388 

$257 

52.2 

27 

I 

$817 

45.2 

70.7 

.417 

$273 

$3523 
-  905 
$2618 

43.9 

II 

$781 

48.0 

74.0 

.370 

$301 

41.7 

36 

I 

$793 

58.1 

75.8 

.408 

$295 

$4573 

-1115 

$3458 

33.8 

II 

$767 

60.4 

78.5 

.360 

$318 

31.6 

50 

I 

$802 

73.4 

81.3 

.394 

$281 

$6205 

-1441 

$4764 

22.0 

II 

$786 

74.8 

83.8 

.346 

$300 

20.7 

75 

I 

$907 

91.2 

88.9 

.373 

$167 

$9120 

-2024 

$7096 

8.2 

II 

$900 

91.6 

90.9 

.327 

$183 

7.5 

100 

I 

$1149 

98.9 

94.9 

.356 

$  o 

$12035 
-  2200 
$  9835 

1.9 

II 

$1148 

98.9 

95.3 

.314 

$  0 

1.8 

Table  B-3.  Summary  of  simulation  results  for  house  size  I,  insulation 
type  I ,  3  cover  collector. 
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AC 

isq 

m) 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

°/o 

SOLAR 

HOT 

WATER 

SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

(million 

kilojoules) 

9 

I 

$586 

17.6 

51.9 

.350 

$98 

39.3 

II 

$563 

m 

54.7 

.300 

$112 

mm 

18 

I 

$540 

37.6 

64.3 

.259 

$142 

$2164 
-  633 
$1531 

29.5 

II 

$529 

39.6 

64.7 

.240 

$147 

28.5 

B 

I 

$540 

53.3 

69.3 

.213 

$140 

$3059 
-  812 
$2247 

22.6 

II 

$534 

54.9 

70.6 

.210 

$142 

22.1 

36 

I 

$563 

65.2 

73.9 

.187 

$116 

$3953 
-  991 
$2962 

mm 

II 

$553 

67.2 

74.3 

.184 

$122 

16.8 

50 

I 

$621 

79.2 

78.7 

.153 

$55 

$5345 

-1269 

$4076 

11.6 

II 

$615 

80.1 

79.3 

.149 

$59 

11.1 

B 

I 

$779 

92.9 

84.7 

.118 

$  0 

S  7830 
-1766 
$6064 

5.5 

II 

$774 

93.3 

86.1 

.092 

$  o 

5.1 

100 

I 

$990 

97.7 

88.6 

.094 

$  o 

$10315 
-  2200 
$  8115 

3.0 

II 

$988 

97.5 

89.8 

.063 

$  0 

2.9 

Table  B-4.  Summary  of  simulation  results  for  house  size  I,  insulation 
type  II,  1  cover  collector. 


AC 

•,sq 

m> 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

(million 

kilojoule^ 

I 

$545 

24.3 

58.9 

.423 

$138 

$1347 
-  404 
$  943 

3  5.6 

II 

$526 

25.6 

61.7 

.  383 

$148 

34.0 

1  Q 

I 

$488 

49.1 

79.3 

.3  54 

$193 

$2319 

24.0 

1  O 

II 

$469 

51.4 

73.5 

.326 

$205 

—  oo4 

$1655 

22.5 

■ 

I 

$476 

68.5 

78.2 

.325 

$202 

$3291 
-  858 
$2433 

15.6 

II 

$458 

70.5 

81.3 

.276 

$214 

14.1 

36 

I 

$494 

82.3 

83.2 

.303 

$181 

$4263 

9.6 

II 

$484 

83.2 

85.5 

.2  54 

$186 

“XU  DD 

$3210 

8.8 

CA 

I 

$569 

94.3 

89.3 

.277 

$100 

$5775 

-1355 

$4420 

4.1 

DU 

II 

$561 

94.7 

91.6 

.227 

$104 

3.5 

B 

I 

$782 

99.9 

96.4 

.243 

$  0 

$8475 

-1895 

$6580 

0.7 

II 

S  780 

99.8 

97.5 

.200 

$0 

0.6 

100 

I 

$1056 

100.0 

99.2 

.224 

$0 

$11175 
-  2200 

S  8975 

0.1 

II 

100.0 

99.3 

.189 

$0 

0.1 

Table  B-5.  Summary  of  simulation  results  for  house  size  I,  insulation 
type  II,  2  cover  collector. 
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AC  HEAT  ANNUAL  % 


$489 


$470 


$4e2 


$466 


S5C3 


$4  94 


$590  97.2 


$586  97.2 


$835  100.0 


S835  100.0 


$1155  100.0 


$1155  100.0 


=>/ 

/o 

SOLAR 

HOT 

WATER 

SOLAR 

SYSTEM 

ANNUAL 

EFF. 

59.3 

.451 

63.1 

.407 

51.3 


53.2 


71.2 


72.9 


85.7 


86.4 


72.6 


5.9 


80.0 


82.7 


85.3 


87.6 


92.6 


94.0 


99.5 


99.5 


99.9 


99.9 


$191 


$203 


$196 


$205 


$170 


5 


$  77 


$78 


$  0 


$  0 


$  0 


$  0 


COST  (mil lion 

kilojoule^ 


35.4 


33.4 


$2474 
-  695 


$1779 


$3523 
-  905 


$2618 


$3458 


$6205 

-1441 


$4764 


$9120 

-2024 


S  7096 


$12035 
-  2200 


$  9835 


22.9 


21 . 3 


14.3 


13.0 


8.0 


7.2 


2.5 


2.1 


0.9 


0.9 


0.2 


0.2 


Table  B-6.  Summary  of  simulation  results  for  house  size  I,  insulation 
tyoe  II,  3  cover  collector. 
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AC 

Csq 

m) 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

S  PACE 
HEAT 

% 

SOLAR 

HOT 

WATER 

SOLAR 
SYS  T  ELM 
ANNUAL 
EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 
ENERGY 
(mill  ion 
kilojoules) 

Q 

I 

$502 

25.2 

55.6 

.342 

S  2  9 

$1270 

26.5 

11 

$480 

26.7 

58.4 

.287 

$36 

-  J81 

$  889 

24.8 

18 

I 

$473 

49.  9 

68.7 

.  246 

$56 

$2164  | 

-  633 
$1531 

18.0 

II 

$456 

53.1  | 

69.4  1 

.220 

$60 

16.6 

■ 

I 

$489 

68.3 

74.7 

.197 

i 

$39 

$3058 
-  812 
$2246 

12.4 

■ 

II 

$475 

71.1 

75.7 

.185 

$40 

11.3 

36 

I 

$528 

80.5 

78.6 

.166 

$  0 

$3953 
-  991 
$2962 

8.7 

II 

$514 

83.1 

30.7 

.153 

$  o 

WMM 

50 

I 

$613 

92.1 

83.7 

.132 

S  0 

$5345 

-1259 

$4076 

5.0 

II 

S  600 

93.3 

86.9 

.098 

$  o 

3.9 

75 

I 

3  815 

98.6 

39.2 

.096 

$  o 

$7830 
-1766 
$  6064 

mm 

II 

$807 

99.2 

91.9 

.053 

$  0 

mm 

100 

I 

$1045 

99.9 

92.2 

.076 

$  o 

$10315 

1.5 

II 

$1038 

100.0 

94.9 

.036 

$  0 

6  JU 

S  8115 

1.0 

Table  B-7.  Summary  of  simulation  results  for  house  size  I,  insulation 
type  ill,  1  cover  collector. 
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AC 

(sq 

m) 

HEAT 

LOSS 

MODEL 

A  ACTUAL 
ENERGY 
COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

SOLAR 

C  ’/  C  ^  T  V 

EF  F 

ANNUAL 

SAYINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 
ENERGY 
(mill  ion 
kilojoule^ 

9 

I 

$465 

33.6 

63.5 

.413 

B 

$1347 

-  404 

$  943 

23.0 

II 

$442 

35.3 

66.9 

.372 

$74 

21.1 

18 

I 

$423 

65.5 

76.4 

.338 

Si  04 

$2319 
-  664 

12.3 

II 

$401 

68.7 

79.8 

.  301 

$112 

$1655 

11.0 

27 

I 

$440 

85.0 

83.8 

.304 

$84 

$3291 
-  858 

$  2433 

6.7 

II 

$423 

87.7 

87.2 

B 

$88 

5.2 

36 

I 

$488 

95.5 

89.6 

.279 

$32 

$4263 

3.1 

II 

$475 

96.8 

93.1 

.223 

B 

-IUjj 

$3210 

2.1 

50 

I 

$604 

99.8 

95.5 

.251 

.-  0 

$5775 

-1355 

$4420 

0.9 

11 

$600 

190.0 

98.0 

.199 

$  0 

mem u 

75 

I 

$847 

100.0 

99.5 

.222 

so 

$8475 

-1395 

$6580 

0.1 

II 

$846 

100.0 

99. 9 

.180 

$0 

0.2 

190 

I 

$1127 

100.0 

99.99 

.  209 

$0 

$11175 
-  2200 

$  8975 

0.005 

II 

100.0 

99.99 

.171 

$0 

0.001 

Table  B-3.  Summary  of  simulation  results  for  house  size  I,  insulation 
type  III,  2  cover  collector. 


>  0 


T; 

aoetual 

ENERGY 

COST 

"  ■'  7  '""1 

SOLA? 
c  p\  r£ 

0/ 

/O 

S  3  LA  R 

n  j  i 

WATER 

SOLAR 

SVS't-v 

ANNUAL 
EF  F 

'  •  *  * 

S  A7 I  MGS 

c  OLA?. 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

Cm  i  1 1  ion 
v  ^  o  c  i  cu  1  ^  si 

i 

S4S6 

34.6 

64.3 

.444 

B 

31424 
-  427 

3  937 

2  2.5 

3 

n 

3440 

B 

68.6 

.400 

Bi 

2  0.4 

■ 

i 

B 

B 

78.1 

.406 

3100 

$2474 
-  695 
31773 

11.3 

■ 

ii 

ESI 

71.2 

81.9 

.355 

3109 

10.0 

I 

ESI 

88.4 

36.2 

.  384 

Bi 

$3523 
-  90  5 

32619 

5 . 4 

/  / 

ii 

$432 

90.3 

99.6 

.334 

373 

4.2 

i 

$503 

97.6 

92.3 

.365 

314 

34573 

-1115 

33453 

1 . 3 

io 

11 

3494 

98.7 

95.4 

.316 

317 

1.2 

53 

1 

$  63  5 

100.0 

33.9 

.344 

S  0 

-.36205 

-1441 

34764 

0.2 

11 

$633 

100.0 

39.3 

.  300 

3  0 

0.05 

r 

3907 

100.0 

99.39 

.324 

3  0 

39120 

-2024 

$7096 

0.002 

7  o 

ii 

3907 

100.0 

39.99 

.  289 

$  0 

0 . 002 

i 

31223 

100.0 

39.99 

.324 

S  0 

312035 
-  2200 

3  9835 

0.001 

ii 

B 

100.0 

99.99 

.  237 

3  0 

0.001 

Table  3-9.  Suiamary  of  simulation  results  for  house  size  I,  insulation 
type  III,  3  cover  collector. 


Table  B-10.  Summary  of  simulation  results  for  house  size  II, 
insulation  type  I,  1  cover  collector. 


*  J 

MOTEL 

ACTUAL 

ENERGY 

COST 

C  r  IT 

/O 

SCLAR 

WATEP 

SOLAR. 

SYSTEM. 

ANNUAL 

ANNUAL 

SAVINGS 

SYSTEM 

AE7<  IZIAPY 

Zl’ZPGY 
-  *  1  '  ~ 

V* 

12 

I 

S1237 

13.2 

52.5 

.427 

mm 

S  1 6  7  1 

■ZB 

II 

S  12  3 1 

B 

BB 

.  392 

jjliSjjj 

$  1 1 7  7 

9  2.63 

24 

I 

29 . 4 

B 

sa 

3326 

3  2  9  6  " 

-  "92 

7  5.52 

II 

31133 

gg 

64 . 1 

.353 

BB 

3  217  4 

75.4 

36 

I 

m 

B 

6  3.2 

.351 

3377 

$  4  2  6  i 

_  -  r.  -  o 

6  4.1 

II 

$1036 

44.3 

70.9 

.325 

3409 

3  2  2  13 

61.1 

43 

T 

31039 

54  . 1 

72.3 

BB 

340  3 

7  5  5  5  9 
- 1 3 12 
34247 

51.4 

II 

31064 

5  5.7 

75.9 

.299 

BB 

4  9.4 

i 

7  5 

I 

31105 

73.7 

80.8 

.  30  9 

3  3  36 

33475 
-139  5 

30  .  3 

II 

310  89 

74.9 

32.7 

.262 

3  5  5  3  0 

29.0 

.2  0 

T 

31219 

36.5 

3  5.5 

.239 

BB 

$  11 17  3 

-  y  0  f  1  n 

11 

31219 

86.5 

36 . 3 

.224 

32  34 

HU 

16 . 6 

Table  3-11.  Summary  of  simulation  results  for  house  size  II 
insulation  tvoe  I,  2  cover  collector. 
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m ) 

nr  ^ 

1.0  s  s 

MODEL 

ANNUAL 

ENERGY 

COST 

n/ 

SOLA?. 

S  PACE 
y?:  — 

3/ 

/  3 

SOLA.? 

SOLA?. 

SYSTEM 

ANNUAL 

rpr 

ANT7UAL 
SA.Y  INGS 

SOLA? 
ST'S  TEM 

auxiliary 
ENERGY 
(mi  1 1  ion 
.<  1 1  o  3  tu  1  e  sj 

p 

12 

T 

X 

$1238 

13 . 3 

52.3 

.457 

mm 

317  7  4 

96 . 6 

II 

$1235 

15 . 3 

5  5.6 

.413 

$2  39 

92  ,  2 

24 

I 

$1134 

30 . 3 

63.9 

.430 

$319 

77.5 

II 

$1137 

32.4 

66.3 

.  338 

$349 

73.7 

36 

J 

$1123 

44 . 1 

70.2 

.418 

$371 

S4573 
- 1 1 1  5 
$  3453 

62 . 2 

X  -L. 

$1039 

46.4 

73.6 

.370 

$  4  0  6 

5  3.9 

0 

m 

m 

56.2 

75.1 

.  409 

$39  7 

$5972 
-1394 
$4  577 

49.0 

ii 

$1067 

58.2 

B 

.  361 

$431 

m 

7  - 

$i!20 

77.1 

82.6 

.390 

m 

$9120 

-2024 

$7096 

2  6.6 

ii 

$1111 

77.7 

84.1 

.343 

$392 

25.3 

10  0 

r 

$1277 

39 . 5 

87.9 

.  375 

$203 

$12035 
-  2200 
$  98  3  5 

13 . 1 

m 

$1276 

39 . 5 

39 . 4 

.  330 

$225 

13.0 

Table  3-12.  Summary  of  simulation  results  for  house  size  IX, 
insulation  type  I,  3  cover  collector. 
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m) 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

S  PACE 

SOLAR 

a.n"1 

WATER 

SOLAR 
SYS  TEM 
ANNUAL 
EFF  . 

SOLA? 

qyC-T  w 

COS  T 

AUXILIARY 

ENERGY 

(mill  i or. 

>  i  lo  j  cu  le  s! 

12 

I 

$794 

17.5 

52.3 

.  344 

$139 

$1563 
-  4  70 
$1093 

m 

II 

$772 

19.0 

54.8 

.292 

$154 

—  ”  ^ 

3  2.6 

24 

I 

$742 

36.7 

63.6 

.255 

$189 

$276  1 

-  752 

$2008 

41.4 

II 

$734 

38.0 

64  .  5 

.242 

$193 

40.3 

36 

I 

$746 

51.1 

68.4 

.215 

$184 

$3953 
-  991 
$2962 

32 . 5 

II 

$738 

53.0 

70 . 1 

.211 

$191 

31.9 

43 

I 

$768 

6  3.5 

12. a 

.  191 

$159 

$5146 

25.2 

II 

$765 

64.6 

B 

.  183 

$164 

-1229 

$3917 

KSB 

75 

I 

$386 

81.7 

79.8 

.  143 

$35 

$7830 

-1766 

$6064 

■SB 

II 

$384 

82.0 

80 . 2 

.140 

$42 

14.1 

100 

I 

$1052 

91.3 

34.1 

.  119 

$  0 

$10315 

8.1 

II 

$  1 0  5 1 

91.6 

34.3 

.  101 

$  0 

-  2  2  00 
$  8115 

3.0 

Table  B-13.  Summary  of  simulation  results  for  house  size  II, 
insulation  type  II,  1  cover  collector. 
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ENERGY 
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/O 

SOLAR 

S  PACE 
HEAT 

a/ 

SO  LAP- 
HOT 
WATER 

SYSTEM 

ANNUAL 

ITTT  C 

Wk  k  • 

ANN~UAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

Cm  1 1 1  l  dp. 

<  1  1  C  '  c  u  1  s  s) 

12 

i 

$744 

23.7 

5  8.3 

.417 

$189 

$1671 

49.6 

ii 

BB 

24.7 

61.5 

.380 

$202 

$1170 

4  3.0 

24 

i 

$571 

47.3 

71.0 

.  353 

$259 

$2967 
-  793 
$2174 

34.0 

ii 

$  6  5  5 

49.2 

73.1 

.  325 

$272 

32 . 7 

36 

i 

$  6  5  5 

66  .  3 

77.6 

.  326 

$271 

$4263 

-1053 

$3210 

22 . 3 

ii 

$638 

67.7 

80 . 4 

.279 

$287 

21.3 

48 

i 

$676 

80.0 

82.4 

.  306 

$246 

$  5  5  5  9 
-1312 
$4247 

14.5 

ii 

$665 

80 . 7 

84.7 

BB 

$253 

13.6 

75 

X 

$324 

96.2 

90.7 

.270 

$89 

$8475 

-1895 

$6580 

4.2 

ii 

$822 

95.8 

92.4 

.223 

$94 

4.0 

100 

i 

$1069 

99 . 7 

95.7 

.246 

$  0 

$11175 
-  ?  ?  n  n 

1.2 

ii 

99 . 5 

96.4 

.  204 

$  0 

$  3975 

1.2 

Table  B-14.  Summary  of  simulation  results  for  house  size  II, 
insulation  type  II,  2  cover  collector. 


ai 


AC 

sc 

-•  1 

HEAT 

LOSS 

MODEL 

A'CrUAL 
EM EH GY 

rpc’7' 

a/ 

/o 

SOLAR 

S  PACE 
HEAT 

SOLAR 

SYSTEM 

ANNUAL 

a:cc:al 

SAYINGS 

SOLA? 

system 

COS  T 

AUXILIARY 

ENERGY 

Cm  i  1 1  icr. 

V. lo]  cu  le  s) 

12 

I 

$750 

23.8 

59.2 

.449 

$182 

$1774 
-  5  2  2 
$1242 

49 . 4 

II 

$727 

H9 

62.8 

.  406 

$199 

47.5 

24 

I 

$677 

49 . 4 

71.9 

.417 

$252 

$3173 

-  ft  1  ft 

32.9 

II 

$653 

51.2 

IS 

.  369 

$272 

$2333 

31.0 

36 

I 

$661 

69.2 

79.4 

.400 

m 

■B 

20.8 

II 

$645 

70 . 3 

82.2 

.  351 

$279 

19.5 

43 

I 

$638 

83.3 

84.4 

.  386 

$233 

$5972 

-1394 

$4578 

12.3 

II 

$678 

34.0 

36.4 

.  339 

$243 

11.5 

75 

I 

$362 

93.3 

94.0 

.  359 

$47 

$9120 

-2024 

$7096 

2.3 

II 

$361 

98.1 

94.8 

.  316 

$51 

2 . 3 

100 

I 

$1158 

100.0 

99 . 2 

.342 

$  0 

$12035 
-  2200 

warn 

II 

$1159 

100.0 

93.9 

.  303 

$  0 

$  9835 

Table  B-15.  Summary  of  simulation  results  for  house  size  II 
insulation  type  II,  3  cover  collector. 
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■-> 

n^.A  * 
LOSS 

y.oz  EL 

r 

-r  •  n  r 

— . . «  —  r'.'  j  i 
q 't 

1 ' 

0/ 

."3 

SOLAR 

W  A  7  E  P. 

S  C  LA  R 

S  Y  S  T  EM 
ANNUA L 

H 

S  0  LA  R 
SYSTEM 
COST 

AUXILIARY 
ENERGY 
(m  1 1 1  ion 
:  1  c  ]  c  c :  e  si 

12 

i 

$679 

24.6 

5  6.0 

.335 

$37 

$1563 
-  470 

36  .  5 

ii 

$653 

25.9 

yryi 

.282 

$46 

$1093 

34.4 

24 

i 

$652 

m 

68.4 

.242 

$63 

$2761 
-  752 
$2009 

25.4 

ii 

$629 

51.0 

68.9 

.221 

$70 

2  3.6 

35 

i 

$669 

66.1 

74.2 

.  195 

m 

17 . 7 

n 

$651 

68.9 

7  5.0 

.  183 

16.2 

48 

i 

$716 

73.3 

78.2 

.  154 

$  0 

$5146 

-1229 

$3917 

12 . 5 

ii 

$701 

80.7 

79.8 

.  156 

$  0 

11.2 

75 

i 

$387 

Hi 

84.9 

.  122 

S  0 

$7830 

-1756 

$6064 

5.8 

ii 

$873 

94.9 

37.7 

.085 

$  0 

4 . 6 

100 

i 

$1099 

98.0 

88.3 

.098 

$  0 

$10315 

3.4 

ii 

S 10  39 

98.6 

91.2 

.056 

$  0 

£•  Cm  vJ  U 

$  8115 

2 . 5 

Table  B-16.  Summary  of  simulation  results 
insulation  type  III,  1  cover  collector. 


for  house  size  II 
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c  ~ 

■  ■  / 

ACTUAL 

ENERGY 

r*  >▼« 

o/ 

S  O  LA  R 
SPACE 
HEAT 

3/ 

SOLAR 

HOT 

WATER 

p 

I 

ANNUA  L 
SAYINGS 

w  ^  T 

S  Y  S  7Z  V. 

A'_-y  -  t  ;  p  -_ 

ENERGY 

...  •  -  -  -  •  o  cv 

12 

i 

S631 

B 

63.5 

.  407 

$34 

$16  71 

J  i  I 

ii 

$603 

33.9 

66.7 

.  369 

$93 

~  D  Ji 
$1170 

30.0 

i 

$532 

63.5 

75.9 

.338 

$130 

$2967 
-  793 
$2174 

13.4 

ii 

uss 

65.9 

79 . 1 

.302 

$139 

16.2 

36 

i 

$600 

83.2 

83.2 

.  306 

$10  7 

$4263 

-1053 

$3210 

9 . 7 

ii 

$530 

86.3 

85.3 

.253 

$112 

3.1 

43 

i 

$660 

94.1 

88.3 

.  282 

$42 

$5559 

-1312 

S4247 

4.7 

ii 

$646 

95.2 

91.9 

.  227 

$42 

3 . 5 

r . 

7  ^ 

7 

$337 

100-0 

96.6 

.244 

$  0 

$8473 

-1393 

S6330 

0.3 

/  0 

II 

$331 

100.0 

9  3.5 

S  0 

0  .  4 

10  0 

I 

$1160 

100.0 

99.3 

.  224 

t  A 

J 

3  1 1 1  ~  5 

: .  2 

II 

3 11 5  9 

100-0 

99.3 

.  182 

$  0 

KbhJ 

j  •  *. 

Table  B-17.  Summary  of  simulation  results  for  house  size  II, 
insulation  type  III,  2  cover  collector. 


rr; 

-  *cc 

y.ossL 

amt.'l 

entity 

1 1 c  T 

■5/ 

£  £  LA? 

S  PACE 

utt  2,  -*• 

SOLA? 

C  **  ^  It 

S  Y  S  T  EM 

2,  ' '  r 

2 VS TIM 

ALT  ILIA? 7 
ENTITY 

Cm  ill:  on 
<  .It;  :  .  le^ 

1 2 

2 

3  6  32 

33.3 

64.  3 

.441 

EH 

3  1 7  7  4 

31.3 

II 

3  6  0  5 

34.9 

63.  3 

.  40 0 

EH 

$1242 

_ 

2  9.0 

24 

I 

3  5  3  6 

6  6.0 

77 . 5 

.  40  6 

3125 

m 

/  .  0 

II 

■H 

63.7 

KB 

.  356 

EH 

'  "7 

^  “t  .  / 

36 

r 

3  610 

36  .  3 

B 

.335 

396 

$4  573 
-1115 

3  3  4  5  3 

3.2 

ii 

35  91 

38.3 

33 . 5 

.336 

$10  0 

6  .  6 

B 

i 

3679 

D 

9  2.0 

.  367 

3  21 

35972 

-1394 

3  4  5  7  3 

HBH 

ii 

in 

97 . 3 

94.5 

.313 

EH 

2 . 1 

/  j 

i 

be 

100.0 

9  9.6 

.340 

$  0 

m 

0 . 9 

r  t 

m 

10  0 . 0 

99.9 

.297 

S  0 

0  .  3 

r 

S 12  5  9 

HE 

Bi 

.323 

3  0 

3120  3  5 
-  22  0  9 

3  933  5 

j  .  Z  2 

' 

i.  i. 

100.0 

99.9  9 

.  290 

3  0 

3 . 3  2 

Table  B-l-3 
insulation 


Summary 
type  III, 


of  simulation  results 
3  cover  collector. 


for  house 


size 


II 


Table  B-19.  Summary  of  simulation  results  for  house  size  III, 


insulation  type  I,  1  cover  collector. 


a>|  vc» 


AD-AQ92  436 


AIR  force  INST  OF  TECH  WRlSHT-PATTERSON  *FB  0H  F/r 

ENVIRONMENTAL  influences  IN  THE  SIMULATION  OF  a  SOLAR  SPACE  HEA--ETC(U) 

UNCLASSIFIED  AFIT-CI-80-20T  NL 


END 
1-8 1 


DTIC 


AC 

(sq 

m) 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

15 

I 

$1031 

16.1 

52.1 

II 

$1004 

18.3 

54.6 

30 

I 

$965 

35.4 

62.8 

II 

$964 

36.1 

64.0 

4 

I 

$967 

49.3 

67.6 

II 

$965 

50.4 

69.6 

fin 

I 

$996 

61.7 

71.6 

u  u 

II 

$988 

61.8 

73.4 

75 

I 

$1045 

70.7 

75.5 

II 

$1042 

71.1 

76.3 

100 

I 

$1174 

82.3 

80.1 

II 

$1170 

82.4 

80.1 

Table  B-22.  Summary  of  simulatior 
insulation  type  II,  1  cover  collec 


: 

i 


SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

(million 

kilojoule^) 

.340 

$169 

$1866 
-  560 
$1306 

71.8 

.290 

$196 

69.6 

.253 

$238 

$3357 
-  871 
$2486 

55.1 

.244 

$235 

55.0 

.219 

$231 

$4848 

43.8 

.210 

$238 

-1170 

$3678 

43.6 

.195 

$207 

$6339 

34.8 

.183 

$209 

-1468 

$4871 

34.0 

.169 

$150 

$7830 

27.2 

.163 

$160 

-1766 

$6064 

26.9 

.139 

$17 

$10315 

18.0 

.138 

$28 

17.8 

results  for  house  size  ill 


AC 

(sq 

m3 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

I 

$963 

22.7 

58.4 

ID 

II 

$947 

23.5 

61.1 

30 

I 

$878 

45.6 

70.3 

II 

$859 

47.2 

72.5 

45 

I 

$850 

64.2 

76.8 

II 

$835 

65.1 

79.5 

60 

I 

$870 

77.8 

81.8 

II 

$865 

77.8 

83.8 

75 

I 

$929 

87.2 

85.6 

II 

$925 

87.4 

87.3 

100 

I 

$1118 

95.7 

91.0 

II 

$1116 

96.4 

92.2 

Table  B-23.  Summary  of  aim' 
insulation  type  II,  2  cover 


!9 


SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

(million 

kilojoules) 

.414 

$237 

Uli 

I 

65.4 

.378 

$253 

$1396 

64.1 

.354 

$318 

$3615 
-  923 
$2692 

46.0 

.326 

$342 

44.4 

.328 

$343 

$5235 

-1247 

$3988 

31.3 

.284 

$365 

30.0 

.308 

$318 

$6855 

-1571 

$5284 

20.4 

.261 

$333 

20.0 

.291 

$256 

$8475 

-1895 

$6580 

12.7 

.244 

$270 

12.5 

.268 

$60 

$11175 
-  2200 
$  8975 

5.2 

.223 

$71 

5.1 

n  results  for  house  size  III, 


ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

$973 

22.8 

58.9 

$949 

24.0 

62.5 

$887 

47.1 

71.1 

$860 

48.9 

74.9 

$859 

66.7 

78.5 

$845 

67.5 

81.4 

$885 

81.0 

83.7 

$880 

81.0 

85.5 

$955 

90.1 

88.1 

$952 

90.4 

89.9 

$1190 

98.0 

94.3 

$1187 

98.6 

94.8 

Table  B-24.  Summary  of  simulatioi 
insulation  type  II,  3  cover  collec 

i 

i 


0 


SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

(million 

kilojoule^) 

.448 

$227 

$2124 
-  625 
$1499 

65.2 

.405 

$251 

63.3 

.418 

$310 

m 

44.7 

.370 

$340 

42.5 

.402 

$333 

$5622 

-1324 

$4298 

28.9 

.353 

$354 

27.8 

.388 

$302 

$7371 

-1674 

$5697 

17.7 

.341 

$317 

17.5 

.376 

$226 

9.9 

.329 

$241 

9.7 

.358 

$  o 

$12035 
-  2200 
$  9835 

2.8 

.315 

$  o 

2.6 

results  for  house  size  III 
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AC 

(sq 

m) 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

Of 
i O 

SOLAR 

HOT 

WATER 

SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

(million 

kilojoule^ 

15 

I 

$869 

24.0 

56.0 

.331 

$44 

$1866 
-  560 
$1306 

47.3 

II 

$841 

25.1 

58.7 

.279 

$55 

45.1 

30 

I 

$837 

46.8 

68.1 

.239 

$74 

$3357 
-  871 
$2486 

33.4 

II 

$815 

49.1 

68.4 

.223 

$80 

31.6 

45 

I 

$858 

64.5 

73.8 

.194 

$50 

$4848 

-1170 

$3678 

23.6 

II 

$842 

66.4 

74.1 

.191 

$52 

22.3 

60 

I 

$918 

76.5 

77.6 

.165 

$  o 

$6339 

-1468 

$4871 

17.3 

II 

$902 

78.2 

78.7 

.159 

m 

15.7 

75 

I 

$996 

85.7 

80.9 

.144 

$  0 

$7840 

-1766 

$6064 

12.0 

II 

$981 

86.7 

82.9 

.127 

$  0 

10.8 

100 

I 

$1175 

94.1 

86.3 

.117 

$  0 

$10315 
-  2200 
$  8115 

6.9 

II 

$1162 

94.8 

87.7 

.081 

$  0 

5.8 

Table  B-25.  Summary  of  simulation  results  for  house  size  III 
insulation  type  III,  1  cover  collector. 
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AC 

(sq 

m) 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

Cmillion 

kilojoule^) 

15 

I 

$811 

31.7 

63.1 

.403 

$101 

$1995 
-  599 
$1396 

41.8 

II 

$800 

33.1 

66.5 

.367 

$115 

39.2 

30 

I 

$753 

61.5 

75.4 

.339 

$155 

$3615 
-  923 
$2692 

24.5 

II 

$724 

63.7 

78.6 

.303 

$169 

22.2 

45 

I 

$770 

81.3 

82.5 

.308 

$133 

$5235 

-1247 

$3988 

13.5 

II 

$748 

82.8 

85.6 

.256 

$141 

11.7 

60 

I 

$840 

92.8 

87.9 

.284 

$58 

$6855 

-1571 

$5284 

6.7 

II 

$824 

93.5 

91.0 

.231 

$60 

mm 

75 

I 

$949 

97.9 

92.3 

.265 

$  0 

$8475 

-1895 

$6580 

3.2 

II 

$937 

98.4 

94.7 

.213 

$  o 

2.2 

100 

I 

$1202 

100.0 

96.9 

.242 

$  o 

$11175 
-  2200 
$  8975 

0.9 

II 

$1197 

100.0 

98.3 

.195 

$  o 

0.5 

Table  B-26.  Summary  of  simulation  results  for  house  size  III 
insulation  type  III,  2  cover  collector. 
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AC 

(sq 

m) 

HEAT 

LOSS 

MODEL 

ANNUAL 

ENERGY 

COST 

% 

SOLAR 

SPACE 

HEAT 

% 

SOLAR 

HOT 

WATER 

SOLAR 

SYSTEM 

ANNUAL 

EFF. 

ANNUAL 

SAVINGS 

SOLAR 

SYSTEM 

COST 

AUXILIARY 

ENERGY 

(million 

kilojoule^) 

15 

I 

$815 

32.7 

63.8 

.440 

$97 

$2124 
-  625 

41.1 

II 

$783 

33.6 

68.2 

.400 

$112 

$1499 

38.4 

30 

I 

$756 

64.2 

77.0 

.407 

$151 

$3873 
-  975 
$2898 

22.8 

II 

$728 

66.2 

80.5 

.358 

$164 

20.5 

45 

I 

$782 

84.4 

84.6 

.387 

$120 

$5622 

-1324 

$4298 

11.4 

II 

$763 

85.7 

87.4 

.338 

$125 

9.9 

60 

I 

$862 

95.5 

91.0 

.369 

$33 

1 

mu 

mm 

II 

$850 

96.2 

93.5 

.320 

$32 

3.5 

75 

I 

$988 

99.5 

95.8 

.354 

m 

Blf 

■ EbceB 

1.5 

II 

$982 

99.5 

97.3 

.310 

$  o 

1.0 

100 

I 

$1293 

100.0 

99.7 

.338 

$  o 

$12035 
-  2200 
$  9835 

0.09 

II 

$1293 

100.0 

99.8 

.297 

$  0 

0.05 

Table  B-27.  Summary  of  simulation  results  for  house  size  III, 
insulation  type  III,  3  cover  collector. 
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//brooks  Job - * 'b rooks *doud'  »elass-a*time=<4*00> 

/♦Jobparm  d=rmtl 

//#  usr=(dmbrook*/bin/sus*l> 

//  exec  forascla*parm.fort=nosource 
//fort.sasin  dd  # 

real  ms  *  nice  * mcpsi  If  svl  »  If  sv2 
INTEGERS  KDATA(8) 
dimension  sum ( 18) * csti (3 *3) 

data  csti<l»l>*csti(l*2)»csti<l*3)/-26.82*0.0»73»07/ 
data  csti<2*l)*csti(2*2)*csti<2*3)/-38.02*0.0*104*01/ 
data  csti (3* 1 ) *csti <3*2) *csti (3 *3) /-49*97*0* 0*  137*09/ 
print  44 
print  44 
print  161 
c 

c  enter  data  values  for  latitude*collector  tilt*  collector 
c  azimuth  -  all  in  decrees, 
c 

ph=35*  50 
betas'SO.S 
3ama=0*0 
c 

c  enter  values  for  max  collector  size*acmx*  no*  of  dlass  covers*n* 
c  desidn  room  temperature*  and  upper  limit  for  service  hot 
c  water  temperature* 
c 

acmx-125* 

n=2 

tr-293« 16 
tr2=291 *48 
thw=333. 16 
c 

c  enter  values  for  collector  cost/so  m*  storage  cost/s a  m* 
c  piPinS  and  equipment  costs*  eauivalerrt  uniform  annual 
c  cost  factor »present  worth  factor*  and  fuel  costs/million  kJ* 
c 

nm2=n-2 

ccost=100* +8* 6*f loat < nm2) 
stcst=8.0 
ce=375. 
euacf =♦ 11746 
euac3=. 10185 
pwf =8 « 932 
cf=* 000006 
c 

c  enter  flad  for  tape  radiaton  values  on  this  data  tape 
c 

c  itape=0  no  direct  normal  values  on  tape* 

c  itape.ne.O  direct  normal  values  are  available* 

c 

itape*0 

c 

c  enter  sise< isize) »desian  air  temp 

c  (tdesn)*  desidn  wind  speed < »vdesn> »  max  and  min  around  temperatures 
c  (tdmaxftsmin) . 
c 
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insul=2 

isise-2 

tdesn=279.2 

vdesn=ll  .5 

tgmax=294.4 

tgifiin=286.46 

t s in  n  =  (  t s m a x + 1 S in  i  n  )/2. 

tSaiiiP=  <  iSmax-tSmin )  /2 . 

print  159 

write<6?  162)  tdesn? vdesn? timin' tSmax 
print  44 
c 

do  900  nac=l?4 

ac=0 ♦ 2*acmx$ f 1 oat  <  nac ) 

cf uel=cf #3c#euac3#20. 

c  call  subtroutine  initial  to  calculate  initial  one  time  constants 
c  for  the  station?  house?  and  collector  used. 

call  initl <ph?beta?sama?sinph?cosph? tanph?sinb?cosb? 
lsind?  cosS?ac?ms?niCPS  ?us?  tsold?strls?tpp?flo?  rho  ?  vdesn? 
lul ?mce?cpc  ?  f  r ?n?exptrm ?  absor  ?  f low?  f p  ?ex  ?  fc? fee ?strl2?tsol2) 
c 

print  44 
print  44 
print  158 

write (6? 40)  ms?f lo?ex?n?ul ?f r 
c 

c  call  subroutine  design  to  compute  station  specific  house 
c  parameters?  areas  of  house  surfaces?  area-u  value  products 
c  for  exposed  surfaces?  design  heat  loss/deg  c?  and  infiltration 
c  loss. 

cal 1  design < ac  ? beta? isize? insul ?  tdesn ?achs  ?  a  ?b ?  c  ?d?e? 
larn?ars?  ahs  ?  a  raw  ?  aw ins  ?  aws  ?  aw inn ?  awn  ?  swine ?  a wen  ?  aud  ?  vdesn  ? 
la  run  ?  a  rus  ?  a  ruaw ?  auc  ?  auwi n  ?  auw  ?  auf 1 ?  auc2  ?  ausum ?  u  rf  ?  adesn  ? 
led  inf  ?  samin  ?uc?uw?ug?uaw  ?ud  ?uf  1  ?adr  ?  fc  ?  aub  ?  tginin?  vohw? 

1 achw  ?  achs  ?  vol ) 
c 

c  print  all  house  parameters, 
print  44 

print  5? isize?insul 
samin=samin#57.28 
hll=2.#aud/adesn#100. 
hl2=auwiri/«desn#100 . 
hl3=eufl /adesn* 100. 
hl4=auw/adesn*100. 
hl5=auc/adesn*100. 
hl6=adinf /adesn#100. 
htl=hll+hl2+hl4+hl3+hl5+hl6 
print  7 

print  6?hll?hl2?hl3?hl5?hl4?hl6?htl 
print  88 

print  6?a?b?c?d?e?samin 
ssmin=sani  in/57. 28 
print  9 

print  21 ?arn?ars?ahs?araw?3wins?aws?awinn?awn?awine?awen 
print  10 

print  21 ?aud?arun?arus?aruaw?auwin?auw?auf 1 ?auc?ausum 
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print  11 

print  21 >adesn>adinf »  schs»  sc 
print  44 
hl=adesn 
IDAY=0 
c 

c  run  the  model  for  the  design  year 
c 

do  64  J=l»18 
64  sum(J)=0* 

DO  1  J=1 t 365 
IDAY=IDAY+1 
DAY=FLOAT  ( IDAY ) 

c  compute  decl i nation*  cos  end  sin  of  declination*. 
ar<3=  <  360 .  *  <  284 .  +day )  /365 .  )  *  .01745329 
delt=(23.45*sin<ard) >*.01745329 
sdelt=sin<delt) 
cdelt=cos(delt ) 
a  rS2=a  rS- .516426 

c  compute  around  temperature  for  this  day 
ts=tamn+tasmp*sin<  ars2> 
c  compute  basement  losses/hr  for  this  day. 
bls=aub#< tr-ts) 
if (bis. le.O)  bls=0. 

if ( i day. at . 150. and. iday. It. 255)  bls=0. 

DO  1  1=1 f 24 

c  read  weather  and  solar  radiation  parameters  for  each  hour. 
READ  (4*4)  ( KDATA ( nv ) » nv=l » 8 ) 

4  FORMAT (8A2) 

IRHT  =KDAT A ( 4 ) 

M3DA=KDATA ( 1 > 

I T  XME=KDAT  A  <  2  > 

TIME=FL0AT ( ITIME) 

ITEMP=KDATA<5) 
ta=f loat ( itemp>/10. 
ta=ta+273. 16 
IDN=NDATA(3> 
rdn=f loat ( idn) 

DN=FL0AT ( IDN ) 

RHT=FL0AT ( IRHT  > 

IC=KDATA<8) 

CC=FL0AT  < IC ) 

IUIND=KDATA(6> 

UIND=FL0AT ( I WIND) 

I  SF'ED=KDATA  <  7  > 

SFEED=( FLOAT (ISPED) )/10. 

c  set  initial  solair  temperatures  to  dry  bulb  temperature. 
tern=ta 
ters=ta 
teww=ta 
teew=ta 
tesw=ta 
tenw=ta 
alpha=0. 
theta=l. 570716 


c 
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c  if  radiation  values  are  greater  than  Or  compute  the  incidence 
c  angle  on  the  collector  and  all  building  surfaces.  If  the 
c  angle  of  incidence  on  ana  surface  exceeds  90  deg*  set  it  =  90. 
c 

if<idn.le.0.or. idn. gt . 9998 )  rdn=0. 
if<irht.le.0.or.irht.gt. 9998 )  rht=0 . 
if  (  i  rht .  le.O  .or .  i  rii  t  .gt .  9998 )  go  to  8 
c 

call  angle ( sineh * cosph* tanph * sinb*eosb r sing t cosg r 
ldaa  r  time  >  theta  r  alpha  r  angws * answe  *ariSww*  anSnr  * 

1 s am in  r angn w r ph r sina  r costh* cdel t  r sdelt ) 
c 

8  continue 
c 

c  call  subroutine  collector  to  compute  the  solar  radiation  incident 
c  on  each  collector  and  the  useful  enerSa  collected*  if  aria*  ba  each, 
c 

call  col 1 ( ac  *  f r  *  ul » tsold*  ta  * tpp* f lo  *cpc  * tco  *  rdn  t rht  *cc r 
li tape* alpha* theta* beta*  rho*au*ac>  rd*  rdu* i rht *  speed  * n * uls * au2 * 
lcosb » sina  *costh*mce  *expt rm»  absor  *  wind  *fcc*tsol2*ult*fp) 
c 

sum ( 1 ) =sum  < 1 ) +ac 
c 

c  if  the  sun  is  up*  calculate  the  solair  temperature  of  each 
c  exposed  building  surface, 
c 

solgn=0. 

if < alpha . le. 0 .2)  go  to  102 

call  sol  ai  r  <  angws  *angwe  *  anSww  *  arignr  *  theta  *  tern  *  ters*  teew* 

1  teww  *  tesw  *  tenw  *  rdn  >  rdu  *  anS.nw  *  speed ) 
c 

call  shg(angwe»  rdn*  rdu* insul *shgn*moda) 
solSn=solSn+shSn*awine 

call  shg ( angws*  rdn  *  rdu*  insul  * shsri*moda) 
solsn=solgn+3wins#shgn 

call  shg (anaww*  rdn*  rdu* insul *shgn*moda) 
salgn-solSn+shSn#awine 

call  shg  <  angnw  *  rdn  *  rdu  > insul *  shsn*moda ) 
so 1 3n=so 1 gn+shSn#awi nn 
102  continue 
c 

c  set  energa  reauirements  for  service  hot  water  for  this  hour, 
c 

ashw=0. 
oshw2=0 . 

if ( i time. gt .5. and .itime.lt.21)  ashw=vahw#< thw-ts ) 
ashw2=ashw 

sum  (  4 ) =sum ( 4 ) +ashw/ac 
c 

c  compute  building  heat  demand  using  design  heat  transmission 
c  coefficient  (computed  in  subroutine  design), 
c 

hls=0. 

if (ta.lt. tr2)  hls=hl*(tr2-ta) 

aload=hls-strls 

if (ctload. It .0. )  oload=0. 


sum<  9 ) =sum ( 9 )  +aload/ac 
sum  (  5 ) =sum  <  5 ) +hl s/ac 
sum  (  6 ) =sum ( 6 ) +st rls/ac 

commute  storage  temperature  and  storage  losses  for  this  hour, 
call  str3e(«C)f  lo»Qu»aload>  tsold>tsnewf  QduniP»uauXfStrls» 

1  oshw .  c»axw  »  vuhw  .  ms  »  »cps»  us  ) 
sum ( 12 ) =sum (12) +naux/dc 
sum(  13 )  =suin  ( 13)+aaxw/ac 

compute  buildina  heat  demand  as  a  function  of  temperature . 
wind.  and  radiation  influences. 

heatl=0. 

call  demand ( speed. wind  »  f c  » a wen » swine  » a  raw  . uw  »u3 » teew.  aws.  a wins . 
ladrrudfteswf  arsrurf  fters>teww»awn»  awinn  » arn. tenw . tern »  auc2  .  act 
ltppf achw»achs» vol  . bis. aub. heat 1 .winls.hli.ta) 

sum ( 7 ) =sum <  7 ) +heat l/ac 
sum(8)=sum(8)+strl2/ac 
ol=hestl-strl2 
if <ol ♦ It .0. )  al=0. 
if (solan. de.al )  soldn=al 
sum ( 10 ) =sum  < 10 ) +solan/ac 
sum  (11)  =suiti  (11)  +hl  i/ac 
l=al-soldn 
if (al . lt.O. )  al=0. 
sum ( 18 ) =sum (18)+al/ac 

compute  storage  temperature  and  storaae  losses  for  the  expanded 
parameter  model. 

call  strde(ac»flo.«u»al»tsol2»  tsnw2 1  c<dmp2 1  aau::2  .  st  r  12 1 
lashw2 » oa>:w2 » vahw » ms  »  ihcps  »  us ) 
sum  ( 14 )  =sum  (14)  +oau>;2/ac 
sum  ( 15 )  =sum  ( 15 )  +aa;;w2/ac 
sum ( 16)=sum( 16 ) twin 1 s/ac 
sum (17)=sum(17)+bls/ac 
sum ( 2 ) =sum ( 2 ) +au/ac 
sum ( 3 ) =sum ( 3 ) +ou2/ac 

1  continue 

compute  annual  efficiency. eauivalent  uniform  annual  costs, 
percentade  solar,  annual  savinds  for  each  system. 

anef l=sum( 2) /sum < 1 ) 
anef 2= sum ( 3 ) /sum  < 1 ) 
colcs=(ccost+stcst)#ac+ce 

if (coles. at. 2000. )  taxsv=600 . +0 . 2* ( colcs-2000 . ) 
if (coles. le.2000. )  ta>:sv=0 . 3*colcs 
if(colcs.ae. 10000.)  ta;:sv=2200 . 
ccfn=colcs-taxsv 

euac 1 =ccf n#eu3cf + ( sum ( 12 )+sum( 13) >*efuel+csti ( isize. insul ) 
euac2=ccfn*euacr+(sum( 14)+sum( 15) >*cf uel+csti ( isize. insul ) 

phwl  =  < 1 . -sum ( 13 ) /sum ( 4 ) )*100. 
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phw2=(  1 . -sum (  15)/sum(4>  )#100. 
pshtl=< 1 . -sum( 12 ) /sum (9 ) )#100. 
psht2=(l.-sum(14)/sum< 18)  )*100. 
e 

fcstl=cfuel*(sum(9)+sum<4) ) 
f cst2=cf ue 1# < sum  < 13 ) +sum <  4 ) > 
if (festi . le .euacl )  ansvl-O* 
if <fcst2.1e«euac2)  ansv2-0. 
if <fcstl .St.euael >  ansvl=f cstl-euacl 
if (fcst2«St*euac2)  ansv2;=fcst2-euac2 
If svi=ansvl*pwf 
If sv2=ansv2#pwf 
c 

c  print  results 

do  333  Jk=l*18 
333  sum<  Jk)=sum<  Jk)#ac 
print  300 
print  100 

write(6f 99)  (sum(kk) »kk*l r 9) 
print  200 

write<6»99)  (sum <kk ) »kk=10» 18) 

print  44 

cful=cf uel/ac 

print  400»ccost»cful 

print  151 >euacl ranef 1 

print  153»phwl »pshtl 

print  152»euac2»arief2 

print  153»phw2f psht2 

print  163» csti ( isize t insul ) 

print  44 

print  155f colcs»taxsv»ccfn 
print  156»ensvl t lfsvl 
print  157»ansv2» lfsv2 
print  44 
c 

REWIND  4 
900  continue 
c 

21  format < 10<lx»fl2»2)) 

44  format < lx »//> lx rl20( '*')»/) 

5  format < lxf 'house  size  ' » i2»5x t  ' insulation  level  /»i2»/> 

6  format < 10<lXffll»2)) 

9  format( lx»/»6xr 'house  surface  areas' r/t 
16x» 'north  roof  south  roof-ac  house  area  attic  walls' r 
1'  s  windows  south  walls  n  windows  north  walls'  r 
1'  end  windows  end  walls  ') 

88  format<  lx»/»5x»  'at  roof  ht  hm-roof  In  crdis  to  pk'» 

1'  drs-roof  In  erhwd-c  slope  n-roof') 

7  format  <  lxf/»6xf '%hl  doors  '/.hi  windows  Zhl  basemt  Zhl  ceil' 
1» 'ns  %hl  walls  %hl  infilt  total  %') 

10  format < lx»/f6x»  'area  u-value  products'  t/t 

16x » '  doors  north  roof  south  roof-ac  attic  walls' »2x» 

1'  windows  walls  floors  ceilinS  sum') 

11  format (lxf/r lx» 'desisri  loss/deS  c  desisfn  infil/c  achs/hr'» 
1'  area  collector') 

99  format < lx» 9(el 1 « 4» lx) */) 
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100  format  <  lx?/?  lx? '  iricid.  rad.'?lx?'ou  coll  l'?3x?'ou  c  oll2'?3x? 

l'hot  water' ?3x? 'Heat  loss  1 '  ?  lx  ?  '  strde  lss  1 '  ?  lxr  'heat  loss  2'?lx? 
l'strse  lss  2 ' ? lx?  'heat  load  1'?/) 

200  format  ( lx?/?  lx?  'solar  Sain  '?  lx? '  inf  ilt  ration '? lx?  'auxil  liars  l'» 
llx?'aux  hw  1 '? lx? 'auxil liars  2'?lx?'aux  hw  2'?4x?'wndow  loss'?2x? 
l'bsmt  loss' ?3x? 'heat  load  2'*/) 

151  formats  lx?/? lx? ' the  standard  system  has  eouivalent  annual  costs' ? 
1'  of'?fl0«2?'  dollars?  with  annual  efficiency  of  '?f8»4?//) 

152  formatdx?/? lx? ' the  expanded  parameter  system  has  eouivalent'? 

1'  annual  costs  of  '?fl0.2?'  dollars?  with  annual  efficiency  of'? 

If 8 . 4  ?// > 

300  format < lx? 'annual  energy  totals?'?/) 

400  format < lx? 'annual  cost  and  efficiency  analyses  for  collector  '?/? 
llx? 'costs  of  '?f8»2?'  dollars  per  so  m  and  fuel  costs  of  '?fl0.8? 

1 'dollars  per  kilo Joule? ' ?//) 

153  format  <  lx?  '  it  provides  '?f8.3?'  ’/.  of  hot  water  and  '? 

If8.3?'  %  of  space  heat'?/) 

155  formats  lx? 'the  system  cost  above  base  is  '?fl0.2? 

1'  dollars  '?/?lx?'a  tax  break  of  '?fl0.2?'  dollars'? 

1'  reduces  this  to  a  final  cost  of  '»fl0.2?'  dollars'?/) 

156  formatdx?'  the  annual  energy  savings  for  system  1  are'? 

If 10*2?'  dollars?  and  the  present  worth  of  the  lifetime' ?/?lx? 

1'  savings?  with  money  valued  at  8 V.  is  ' ?f 10.2? 'dollars' ?/) 

157  formatdx?'  the  annual  energy  savings  for  system  2  are'? 

Ifl0.2?'  dollars?  and  the  present  worth  of  the  1 if et ime ' ?/ ? lx ? 

1'  savings?  with  money  valued  at  8%  is  ' ?fl0*2? 'dollars' ?/) 

158  formatdx? 'collector  and  storage  parameters?'?/? 

llx? 'mass  of  storage '? lx? 'mass  flow  rate/sa  m  of  col*'?lx? 
l'heat  exchanger  coef . ' ? lx? 'no*  of  glass  covers'? lx? 
l'col.  u  loss '? lx? 'heat  removal  factor'?/) 

40  format  dx?fl5«l?  6x  ?fl5.2?8x?fl5.4?4x?i8?8x?fll.2?3x?fl5.4?/) 

159  formatdx? 'design  weather  parameters? '  ?/? 
llx? 'design  air  temp' ? lx? 'design  wind  sp'rlx? 
l'winter  ground  temp' ? lx? 'summer  ground  temp'?/) 

161  format (lx? 'Oklahoma  city'?//) 

162  format dx?fl5*2?lx?fl4.2?3x?fl5»2?3x?fl5*2?/) 

163  format (lx? ' the  eouivalent  uniform  annual  cost  of  insulation  is'? 

If 8 . 2 ? '  dollars'?//) 

STOP 

end 

//so.ft04f001  dd  dsn^d. brooks. staa?disp=shr 
// 

/*eof 
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subroutine  design  ( ac  * beta *  d  sice  *  insul  ►  tdesn  *achs  *a*b*c»d*e» 
l3rn»3rs*3hs*ar3w*3wiits*3ws>  aw  inn  *  awn  *  aw  i  lie  *  a  wen  *  aud  *  vdesn  * 

1  a  run  *  a  rus  *  a  ruau * auc * auwiri  *  auw  *  auf  1  *  auc2 *  ausum * u  r  f  *  adesn  * 
lad  inf  *  samintuc  »  uw  »uS  *  uaw  * ud  Mjfl»  adr  *  f  c*  aub * tdmiri * v«hw » 
lschw* achs  * vol ) 
c 


c 

c 

c  / 

c  /  ( 

c  /  tilt 

c  /  /! 

c  /  /  ! 

c  /  /  ! 

c  /  /  ! 

c  /  /  / 

c  /  /  / 

c  /  /hlen 

c  !  I  / 

c  ht !  !  / 

c  !  !  / 

c  ! _ !  / 

c  hwd 

^  ^  ^  ^  y\f  ^  ^  ^  ^  ^  ^  ^  ^  ^  vb  ^  ^  tb  ^  tb  ^  ^  ^  ^  d/  lb  ‘b  ^  ^  ^  ^  ^  ^  «b  ^  df  ^  >b  ^  d<  ^  ^  ^  >b  b  -b  d>  (b  >b  ^b  b  >b  df  «l»  tb  tb  b  ^  >b  ^ 

b  *  *  ^  ^  ^  ^  ^  ^  ^  ^  *p  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^ 

C 

c  this  subroutine  calculates  the  design  heat  loss  for  the 
c  structure  under  consideration  as  a  function  of  the 

c  llector  sizer  the  tilt  of  the  collector)-  the  size 

c  and  dearee  of  insulation  of  the  various  components  of 
c  the  structure*  and  the  design  temperature  and  number  of 
c  airchandes  per  hour  for  the  structure.  The  three  house 
c  sices  (isice  =1*2*  or  3)  are  92.9  sa  in*  134.8  sa  m*  and 
c  ISO  sa  m.  The  three  levels  of  insulation  (insul  =1* 
c  2*  or  3)  are  poor*  moderate*  and  heavy  -  relative  terms, 

c  the  subroutine  has  a  table  of  values  for  u  for  all  the 

c  components  for  each  level  of  insulation*  as  well  as  the 
c  dimensions  of  each  of  the  basic  huuse  sices.  The  roof 
c  dimensions  vary  with  location*  since  the  collector  is 
c  taken  as  an  integral  part  of  the  roof.  These  dimensions 
c  are  calculated  in  this  subroutine. 


c 

c  the  outputs  are  the  areas  of  the  surfaces  of  the  structure* 
c  the  dimensions  of  the  roof  structure  (calculated  here  so 
c  that  the  collector  is  an  integral  part  of  the  roof  at  a 
c  tilt  beta)*  the  design  area-uvalue  products  for  all  of  the 
c  house  components*  and  the  design  infiltration  and  overall 
c  heat  transmission  coefficient  for  the  house, 
c 

«b  lb  kb  kb  d.  kb  kb  ^  ^  ^  ^  ^  lb  ^  ^  ^  kb  lb  a.  ^  ^  ^  ^  ^  ^  lb  |b  lb  b  b  b  da  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b  b 

Q  ^  ^  9  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  *  *  'r  *  *  *  *  *  *  ^  ^  ^  •  T  *  ^  ^  ^  ^  *  *  *  *  ^  ^  ^  ^  ^  ^  4*  ^  ^  ^ 

c  dimension  design  ail — change  values  and  u-values. 
c 

dimension  acw(3) *acs(3) 

dimension  uct(3) *uwt(3) *uf It (3) *ugt(3) *uawt(3) *ur ft (3) *udt(3> 
c 

c  dimension  array  for  house  component  sizes. 
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dimension  sh<3)  , hlen<3)  »hwd<3>  »ht(3)  »f  1  (3)  »awin(3>  , ad(3> 
c 

c  enter  design  u-values  and  aii — change  values  in  table? 
c 

data  uct<l>  ,uet<2) ,uet < 3 > /3 . 883 » 1 . 627 , 0.850/ 
data  uwt(l) ,uwt(2> , uwt < 3 > /5 .070 ,2 . 730 » 1 . 463/ 
data  uf 1 1 ( 1 ) » uf 1 1 < 2 ) , u f 1 1 ( 3 ) /7 ♦ 162 » 2 . 773 ,2.1 68/ 
data  uSt <  1  > » uSt < 2 ) , uSt < 3 ) /25 .409,10.631,9.751/ 
data  uawt ( 1 ) , uawt  <2> , uawt < 3 )/8 . 856 ,8.856,5. 928/ 
data  urft<l> ,urft (2) ,urft<3>/12.338, 12.388,9.975/ 
data  udt <  1 ) » udt ( 2 ) ,udt <3 >/7.563, 6.541 ,5.724/ 
data  sew ( 1 ) , acw<2) ,acw<3)/.030, .0225, .015/ 
data  acs(l> ,acs<2> ,acs<3)/.0216» .0162, .0108/ 
c 

c  enter  dimension  of  each  house  type  in  table? 
c 

data  ah< 1 ) ,hlen( 1 ) ,  hwd< 1) ,ht<l) ,  f  1  ( 1  > , awin ( 1 > , ad ( 1 )/ 
192.9,12.2,7.62,2.44,39.6,9.67,5.0/ 
data  ah (2)  ,hler.<2)  ,hwd<2)  ,ht(2)  ,fl<2)  , awin (2)  ,ad<2)/ 
1134.8,14.25,9.5,2.74,47.5,13.5,6.7/ 
data  ah(3>  ,hleri<3)  »hwd<3> ,ht<3>  ,fl<3)  ,awin<3)  ,ad(3)/ 
1180. ,15.0,12.0,3.05,54.0,17.67,7.31/ 
c 

cpa=l. 20719 
radcon= .01732945 
ars=beta*radcon 
ars2=(90.-beta>#radcan 
snb=sin<arg) 
snb2=sin(ars2) 
c 
c 

i=isise 

in=insul 

d=ac/(hlen( i )*0.8) 

a=snb*d 

c=srib2*d 

e=hwd  < i ) -c 

tarsf=a/e 

samiri=atan  ( tarS) 

b- <  a#a+e#e ) ##0 . 5 

c  compute  area  of  north  and  south  facing  roof  surfaces. 
arn=b#hlen< i ) 
ars=d#hlen(i )-ac 
c  compute  area  of  attic  walls. 
ahs=ah< i ) 
araw=0.5*a*hwd( i ) 

c  compute  area  of  south  facing  windows. 
awins-awin< i >*0 .6 

c  compute  area  of  south  facing  wall  surface. 

aws»hlen( i )#ht< i )-ad( i )/2.-awins 
c  compute  area  of  north  facing  windows  and  walls. 
awinn=awin( i ) * .2 

awn=hlen< i )#ht ( i )-ad< i )/2 . -aw inn 
c  compute  area  of  window  glass  and  walls  on  ends  of  house. 
awine=awin< i ) *♦ 1 
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awerishwd(  i  >*ht<  i  >-auine 
adr=ad( i )/2. 

c  compute  effective  volume-density  of  inside  of  house  for  inf i Itration. 
vol=ahs*ht ( i )#cpa 

c  set  design  infiltration  load  due  to  wind  and  stack  effects. 

3chw-acu( in) 
achs=acs ( in ) 

acha= • 10+achs*  <  291 . 5-tdesn )  +achw*vdesn 
adinf =vol*achg 

c  Set  design  u-values  for  each  surface  for  this  house. 
uc*uct(in> 
uw=uwt (in) 
ug=ugt < in) 
uaw-uawt ( in) 
urf=urft < in) 
ud=udt< in) 
uf l=uflt< in) 

c  compute  design  area-u  value  products  for  this  house. 
aud=adr#ud 
arun=arn#urf 
arus=ars#urf 
aruaw=araw#u3w 
auc=3hs#ue 
auwiri=3win<  i )  *ua 
auwn=awn*uw 
auws=aws*uw 
auwen=aweri*uw 
auw=suwri+3uws+2 .  #auwen 
aub=uf l#f 1 ( i ) 

auf l«aub# ( 291 .5-tamin>/<291 . 5-tdesn ) 
auc2=auc#295. 16 

ausum=auc+arun+arus+2.*3ruaw+3C#urf 
c  calculate  design  attic  temperature 

desnat=ta t ic ( auc2 » tdesn » a  run » tdesn » a  rus r  ac » u  r f  * tdesn  r 
1 tdesn » tdesn  >  a  ruaw » a  ruau » ausum ) 
c  calculate  design  heat  transmission  coefficient. 
adesn=auw+auwin+auf l+auc+2 .  *eud+adirif 
uw=l ./uw-fc 
uaw=l ./uaw-fc 
ud=l ./ud-fc 
ug=l ./ug-fc 
urf =1 ./urf-fc 

c  set  hourly  energy  reouirements  for  hot  water  based  on 
c  house  size, 
c  75kg/person-daw 
c  4  persons  for  size  1 
c  5  persons  for  size  2 
c  6  persons  for  size  3 
kv*isize-l 

vohw*<  <300.+float<kv)*75. >*4.19)/15. 

return 

end 
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function  tat i c ( auc2 » te  rn » a  run  t tens » a  rus » ac » ur f , tpe  ? teew » teww » 
laruae  raruau » ausum) 

c  this  subroutine  calculates  the  temperature  of  attic 
c  usins  the  proceedure  aiven  in  ashraerch  25»p440. 

tatic»<auc2+tern*arun+ters*arus+ac#urf*ters+teew*aruae 

1  Uewwia  ruaw  ) /ausum 
re  turn 
end 
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^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ijj  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  \|r  ^  ^  j|^  ^  ^  jj^  ^  «|^  j^>  ^  ^  j|^  ^  j|/  jjj  ^  ^  ^  j|^  ^  ^  ^  ^  ^  ^|/  <|^  ^  jj^  j||  j|j  ^  j|| 

C 

subroutine  ini  tl  < ph » bets * same » sinph*eosph*  tanph*sinb*cosb* 

1  sins * coss >aoitSf  ihcps » us » tsold » st rls *  Lpp » f lo *  rho r  vdesn * 
lul*mce*cpc*f r*n»exptrm*absor*f low*fp*ex*fc*fce*s trl2*tsol2> 
c 

c  This  subroutine  takes  statiori  latitude*Ph * cul  1  lector  tilt*  c 
c  beta*  and  the  collector  azimuth  angle*  Seme  end  computes  e 

c  the  trigonometric  one  time  calculat ions  needed  for  sub-  c 

c  rout ine *  angle  *  and  initializes  values  for  subroutine  storage* 
c  collector*  and  design. 

c  input  values  for  all  angles  are  in  degrees* 
c 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  «^»  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  v^»  ^  ^  ^  ^  ^  ^  ^  ^  j|^  ^  ^  ^ 

C 

real  ms*mcps*mce 

dimension  f(3)*ult<3) *tmce<3>  *f rt<3> 
data  f <1> *f <2> *f <3>/.95» .98* .99/ 
data  frt<l) *frt<2>  »frt<3)/.786* .824* .901/ 
data  ult(l) *ult<2) *ult (3) /28 . 8 *  14. 4 * 7. 2/ 
data  tmce< 1 ) * tmce<2) * tmce< 3 > /6 .6* 10 . 66*15.22/ 
vcdesn*5« 
radcv=. 01 745329 
epsl=. 1111 
p=ph*radcv 
b=beta#radcv 
g=gama*radcv 
sinph=sin<p) 
eosph=cos(p) 
tanph=tari<p) 
sinb=sin(b) 
cosg=cos<g) 
cosb=cos<b) 
sins=sin(s) 
c  set  storage  mass 
ms-61 .  #ac 

c  set  heat  capacity  of  storage  mass*  deg  c/kJ. 
mcps=l . /(ms*4.19) 

c  compute  area  of  cylindrical  storage  vessel  for  this  storage  mass. 

area=9 . 4247779*  < ( ms/9424 . 7 ) ** . 6666 7) 
c  compute  heat  transmission  coefficient  for  storage  vessel. 
us=l «44*area 

c  initialize  storage  temp*storage  loss*plate  temp*rho. 
tsold-333.16 
tsol2=333. 16 
strls*0. 
strl2»0. 
tpp*273. 16 
rho*0.2 

c  specific  heat  of  water  -  kJ/(kg  -  deg  c>» 
cpc*4. 19 

c  get  collector  heat  loss  coefficient  for  n  covers*  kJ/deS  c. 
ul=ult <n) 

c  get  effective  collector  heat  capacity  -  deg  c/(kJ  -  so  m>. 
mce=tmce<n) 

c  compute  exponential  term  for  plate  temperature  calculations. 
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ara®-<ul/nice> 
expt  rm=eKP <  a  re  > 

c  compute  fr»  heat  removal  factor  for  this  collector* 
fp=f <n> 
f ro=f rt<n) 
eps=.001 
ict-0 

13  continue 

et rm=-2 ♦ * < f p-f ro)/fp 
ict=ict+l 

f rn=<fp#fp/(2**(fp-f ro> ) )#< 1 .-exp(etrro) > 
test=f rn-f ro 
f ro=f rn 

test=abs< test ) 

if < ict .2e. 100)  So  to  14 

if < test.2t.eps)  do  to  13 

14  continue 
f r=f rn 

c  set  mass  flow  rate 

f  lo3Sfp*fp#ul/<2.#cpc#(fp-f  r)  ) 
f low=f lo*ac 

c  calculate  design  outside  film  coefficient 
vc=0 . 2S#vdesn* . 4344881 
fc=67.049*(vc**.605)  +14.0 
fc=l./fc 

c  set  absorbtivits  of  collector  plate  (copper*  tube-in-sheet rblack 

c  non-selective  coating). 
absor=0.95 
return 
end 
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subroutine  angle (  sinph*  cosph  » tanph  t  si  rib  *cosb » sing » cosg» 

Idas  r  time » the  ta  t  alpha  r  andws  *  angwe  t  an^uu  »  ansfnrr 
1  san in  r  ansnu » ph  ?  s me * cost! 1 1  cde  1 1  *  sde  lt> 
c 

C*************************************  *************  t****X*M%XX%MX%X** 
c  This  subroutine  computes  angles  used  in  solar  calculations, 
c 

c  It  reauires  that  the  subr  initial  compute  the  following? 
c 

c  sinph*cosph»tanph  -  the  siriFCos»tan  of  the  latitude, 

c 

c  sins  and  cosg  -  the  sin  and  cos  of  the  azimuth  angle 

c  of  the  surface  for  which  the  incident  solar  angle  is 

c  desired.  Angles  east  of  south  are  positive  and  west 

c  of  south  are  negative, 

c 

c  sinb  and  cosb  -  the  sin  and  cos  of  the  tilt  from  hori- 

c  zontal  of  the  surface  for  which  the  angle  of  incidence 

c  is  desired, 

c 

c  das  and  time  -  daw  is  the  number  of  the  das  of  the  searr 

c  and  time  is  the  time  of  das  on  a  24  hour  clock  basis, 

c 

c  It  returns  the  following? 

c 

c  Theta  -  the  angle  of  incidence  on  the  surface. 

Alpha  -  the  altitude  of  the  sun. 

c 

c 

c  angws  -  angle  of  incidence  on  the 

c  south  wall  of  the  structure, 

c 

c  angww.anawe  -  the  angles  of  incidence 

c  on  the  east  and  west  walls  of  the  structure, 

c 

c  ansnw.ansnr  -  the  angle  of  incidence  on  the  north 

c  wall  and  roof, 

c 

c  these  angles  used  in  later  calculations  of  radiation 

c  and  solair  temperature  and  are  set  euual  to  pi/2  if 

c  >Pi/2.  This  is  done  so  that  calculations  taking  the 

c  cosine  of  the  angle  of  incidence  will  be  0  if  the 

c  angle  is  >  pi/2, 

c 
c 

c ********************  **************** ********************************** 

c 

radcon* . 0 1 745329 
c  compute  hour  angle. 

h*-6.28319*( time-12. )/24. 

coshe»cos(h> 

sinha«sin<h) 

c  compute  cos  of  incidence  angle  on  collector  and  sin  of  solar  altitude. 
costh*sdel t*sinph*cosb-sdel t*cosph*sinb*cosg 

costh“costh+cdelt*cosPh*cosb*cosha+cdel L*sinph*sinb*cosa*cosha 
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costh=costh+cdel  t*sinb*sirid*sinha 
sin3=sinph#sdel  t+cosph<*cdel  ttfcosha 
alpha=arsin(sina) 
if (alpha. le.O. )  do  to  6 

c  compute  incidence  andle  on  all  house  surfaces* 
theta=a  rcos ( costh ) 
an2=ph#radcon+l .5707963 
an=Ph*radeon+samin 
can=cos<an) 
san=siri(an> 
can2=cos(an2) 
san2=sin(an2> 

cosnr=can#cdel  tXtcosha+santfsdel  t 

cosriw=cari2*cdelt*eosha+san2#sdel  t 

cosws=-sdel  t#cosPh+cdelt*s:i.nph>i<cosha 

coswe=cdelt#sinha 

cosww=-cdelt#siriha 

anaww=l .570796 

andue=l .570796 

andws=l .570796 

andnr=l. 570796 

andnw=l .570796 

testw=abs ( cosww ) 

teste=abs  <  cosue ) 

test s=abs ( cosws ) 

testn=abs  ( cosrir ) 

test5=abs  <  cosnw ) 

if < testw.dt . 1 . )  do  to  1 

andwu=a rcos ( cosww ) 

1  continue 

if (teste. dt.l. )  do  to  2 
aridwe=arcos(coswe) 

2  continue 

if (tests. dt.l . )  do  to  3 
andws=arcos ( cosws ) 

3  continue 

if ( testn .dt . 1 . )  do  to  4 
andnr=a rcos (cosnr) 

4  continue 

if (test5.dt« 1 • )  do  to  5 
andnw-arcos ( cosnw ) 

5  continue 

if (cosww. It. 0. )  andww=l .570716 
if (coswe. lt.0. )  andwe=l ♦ 570716 
if (cosws. It. 0. )  andws=l .570716 
if  (cosnw.  It  .0. )  ananw=l .570716 
if (cosnr . It .0. )  andnr=l .570716 
do  to  7 

6  continue 
theta*l. 570716 
andws=l .570716 
andwesl .570716 
andww-1 .570716 
andnwal .570716 
andnr*l .570716 
alpha»0. 
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c 

c 

c 

c 
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sub  routine  col  1  ( ac  ,  f  r  ,  ul  ,  Isold ,  ta  ,  Lpp  fflorcpotcor  rdn  rrhtrc 
1  i  type » alpha  ,  theta  , be  ta ,  ri  10 »  au  » uc  ,  rd  »  rdu , i rh l » speed  ,  n » u Is » au 


lcosb » sina » costh » race  t  e::pt rm f  auso r  » wind »  fee  » tsol2,ult  »fp> 
reel  mce 

dimension  taul <3) ,b(3),slp(3) 
data  taul(l) » taul (2) * taul < 3 ) / . 92 , .86. .75/ 
data  b(l),b(2) > b ( 3 ) /2 . 3657 , 2 . 21 1 43 , 1 . 92857/ 
data  s1p(1) r sip (2 ) , sip ( 3 ) / . 0262S , .02457, .021428/ 

This  subroutine  computes  the  useful  energy  collected  by  a 
water  heating  collector.  The  plate  temperature  must  rise»based 
on  its  heat  capacity , heat  loss  characteristics ,  and  the  incident 
radiation,  at  least  10  dee.  c  above  the  temperature  of  the 
storage  fluid  before  useful  energy  is  collected. 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


the  input  arguments  are, 

1.  ac  -  the  collector  area  in  meters. 

2.  fr  -  the  heat  removal  factor  for  the  collector. 

3.  ul  -  the  heat  loss  coefficient  for  the  collector  in 
kJ/s a  m-c. 

4.  tsold  -  the  previous  hour's  storage  temperature. 

5.  ta  -  the  current  hour's  ambient  air  temperature. 

6.  tpp  -  the  previous  hour'  mean  Plate  temperature. 

7.  flo  -  the  fluid  mass  flow  rate  (kg/hr). 

8.  rdn  -  the  current  value  of  the  direct  normal  component 
of  radiation  ( k J/so  m ) . 

9.  rht  -  the  current  value  of  the  global  horizontal 
radiation  CkJ/sa  m). 

10.  cc  -  the  cloud  cover  in  tenths  of  sky  coverage. 

11.  itype  -  a  flag  denoting  whether  or  not  direct  normal 
radiation  values  are  available, 

itype=0  -  no  direct  normal  values  available. 

direct  component  calculated  from 
global  horizontal  and  cloud  cover, 
itype. ne.O  -  direct  normal  values  are  available 
on  tape  and  are  not  calculated. 

12.  theta  -  the  angle  of  incidence  of  direct  component  of 

radiation,  from  subroutine  angle,  (radians) 

13.  alpha  -  the  solar  altitude  in  radians. 

14.  beta  -  the  tilt  of  the  collector  from  the  horizontal  (rad.). 

15.  rho  -  the  reflectivity  of  the  surface. (dimensionless) 

the  output  arguments  are* 

1.  ou,ctu2  -  the  useful  energy  collected  (kJ/sa  m). 

2.  uls»ult  -  the  collector  losses 

3.  tpp  the  collector  plate  temperature. 

radcon=0. 01745329 

if  (  i  rht  .sit  .0  .and*  irht .  It  .9999)  go  to  8 
rd=0 . 
rdu*0 ♦ 
rdn=0. 
rht=0. 
go  to  6 
8  continue 

if ( itype. ne.O)  go  to  1 
7  continue 
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ckt=0.7*(l.-0.1*cc) 
dbhb=l  .-1  .3575*ekt 
rdu=dbhb#rht 

if < alpha . le . 0 . 2)  So  to  2 
rdn=(  rht-rdu)/sina 
rd=costh#  rdn 
do  to  2 

1  rd- rdn*cos  Lbi 
test2=rdn*sina 

if  ( test2.de.  rht)  do  to  7 
rdu=rht-rdn*sina 

2  continue 

if (alphabet. 0.2)  do  to  6 

rd=0. 

rdn=0 . 

rdu=rht 

6  oc=rd+rdu*0 .5*  ( 1 .+cosb )  +  rht#0 . 5# ( 1 . -cosb ) * rho 
tplate=ta+ (  «c/ul )  -  ( oc/ul-  ( tpp-ta  )  >  #e;;ptnti 
tpp=tplate 

if<oc.le.0«>  do  to  3 
ttest=tplate-10 . -tsold 
if(ttest)  3»4*4 

3  au=0. 
au2=0. 
ul t=0 . 
uls=0. 
do  to  5 

4  continue 
and-theta/radcon 

if < and. le.6G . )  tau=taul<n) 

if ( and . dt .  60 . and. and . It . 90 )  tau=b <n ) -sip  Cn ) *and 

if (and . de . 90 )  tau=0. 

ult=ul# ( tsold-ta) 

auaac*f r# ( ac#tsu#absar-ul t ) 

if(cau.le.O)  ou=0. 

ulsssul#  ( tsol2-ta ) 

au2=ac#f  r*  ( e»c*tau*absor-ul s ) 

if ( qu2 . 1 t . 0 . )  qu2=0. 

5  continue 
return 
end 
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subroutine  strae(ac?  fio  ?  ou  ? a  load  ?  Iso  I  ci  ?  tsnew  ?  oduiup  ?  aaux  ? 

1st  r  1  S  ?  Qshw  ?  Q3MU  ?  VC*hw  ?  ItiS  ?  IHCPS  ?  US  ) 

real  ms?mcps 

c  This  subroutine  computes  the  temperature  of  the  stora.se  fluid  c 
c  for  3  solar  space  heatina  system  based  on  the  useful  enerSs  c 
c  collected  by  a  flat  Plate  collector?  the  heat  loss  character'-  c 
c  istics  of  the  storage  vessel?  the  mass  of  the  storage  fluid?  c 
c  and  a  control  strategy.  Based  on  the  system  characteristics  u 
c  and  control  strategy?  it  also  computes  the  uncontrolled  heat  e 
c  losses  from  storage  to  the  heated  structure?  the  amount  of  e 

c  eneray  which  the  collector  must  ’dump"  if  the  maximum  storage  c 
c  temperature  is  reached?  and  the  amount  of  auxiliary  enersy  c 

c  reauired  if  the  heat  load  cannot  be  met  from  the  solar  system,  c 
data  tmax  ?  tmin/373 . 16  ?302. 16/ 
data  strtp ? thw/2?8 . 16 ? 333 . 16/ 

tsnew=tsol d+mcps*  <  au-us* ( tsold~st rtp ) -eload-oshw ) 

i f  (  tsnew . le . tmax )  tdump=0 . 

if  ( tsnew . at ♦ tmax )  tduc.p- tsnew -tmax 

if ( tsnew .at . tmax )  tsnew=tmax 

if (tsnew . le . tmin )  aaux=«load 

lf(tsnew.le.tmin)  oaxw-oshw 

i f < tsnew . St ♦ tmin )  aaxw=0. 

if (tsnew. St. tmin)  oaux-O. 

i f ( tsnew .It. thw . and . tsnew . at . tmin )  «axw=vehw# ( thw-tsnew ) 
i f ( oshw . le . 0 . )  Q3xw=0. 
if ( aaxw . at . 0 . )  tsnew- tsneu+mcps*aaxw 
strls=us#itsold-strtp) 

if (tsnew. le. tmin)  tsnew=tsold-mcps*s trls+mcps*au 

o dump =t dump /mcps 

au=ou-adump 

tsold=tsnew 

return 

end 
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subroutine  sol 3i  r  ( answs  ?  anSwe  ? andww  ?  ansnr  ?  theta  ?  te  rn? 
Iters  ? teew ? teww ? tesw ? tenw  ? rdn  ? rdu  ? andnw  ? v ) 

this  subroutine  calculates  the  solair  temperature  for 
c  all  of  the  exposed  surfaces  of  the  house  from  the  hourly 
c  radiation  values?  the  snsle  of  incidence  on  each  surface 
c  (calculated  i n  sub r  ansle?  ?  c h e  at* so ro t i v 1 t s  lor  sola?* 
c  radiation  for  direct  and  diffuse  radiation?  and  the 
c  current  value  of  ambient  air  temperature.  If  no  solar 
c  radiation  is  present?  the  solair  temperature  is  the 
c  ambient  air  temperature, 
c 

c  set  combined  radiation  -  convection  heat  transfer  coef. 
c  see  ASHRAE  suide?  chapter  23. 
fcr=42.+15.6*v 

c  determine  product  of  fcr  arid  direct  normal  radiation, 
c  a  value  of  0.7  is  used  for  absorbitvity . 
adrdn=0«7*rdn/fcr 

c  determine  product  of  fcr  and  diffuse  radiation. 
rhdu=rdu*0.7/fcr 

c  determine  solair  temperatures  for  each  surface. 
tern=tern+rhdu+adrdn#cos ( ansnr ) 
ters=ters+adrdn#cos ( theta )+rhdu 
teew= teew+ad  rdn*eos ( anSwe )  +  rhdu 
teww= teww+adrdn#cos  <  andww )  +  rhdu 
tesw=tesw+ad  rdnfccos ( answs )  +  rhdu 
tenw=tenw+rhdu+cos ( anSnw ) #adrdn 
return 
end 


115 


subroutine  demand (vrdir»  f ci a wen  t swine  *  a raw » uw  r  u3 > teew r  aws » aw ins » 
ladr  ?  ud » tesw  t  a  rs  >  u  rf t  Lersi  teww » awri  >  aw  inn  >  arn  >  tenw  r  tern  f  auc2 » ac » 
ltF-p  »  schwi  achs»  vol  »bls  »aub*heatl  r  winls » hi  i  t  ta ) 
t r=293 .16 

c  This  subroutine  calculates  the  building  heat  loss  for 
c  the  house  and  insulation  level  specified  as  a  function 
c  of  temperature*  wind »  arid  radiation  influences* 
c 

c  east  side 

call  fco(v»dir»90. »foc) 
fcoi=l ./foe 
uwve=l ./(uw+fcoi ) 
uawve=l ./(uaw+fcoi ) 
usve=l . /(uS+fcoi ) 
hle=uwve*awen#  ( tr-teew ) 
if < teew.de.tr)  hle-O. 
c  south  side 

call  f co ( v.di r . 180. »foc) 
fcoi=l ./foe 
uwvs=l . /(uw+fcoi > 
usvs=l ./(uS+fcoi ) 
udvs=l ♦ / ( ud+f coi > 
urfvs=l ./(urf+fcoi ) 

hlso=( uwvs*aws+udvs#adr )#( tr-tesw) 
if (tesw.ge.tr)  hlso=0. 
c  west  side 

call  fco(v»dir »27G« »foc) 
fcoi=l ./foe 
uwvw=l ./( uw+fcoi ) 
uawvw=l * /( uaw+f coi ) 
uSvw=l . /(uS+fcoi ) 
hlw=uwvw#awen#( tr-teww) 
if (teww.Se.tr)  hlw=0. 
c  north  side 

call  fco(v»dir.360. »foc> 
fcoi=l./foc 
uwvn=l ./(uw+fcoi ) 
usvn=l ./(uS+fcoi ) 
udvn=l ./(ud+f coi ) 
urfvn=l ./(urf+fcoi ) 
hln=(uwvn*awn+udvn*adr )#( tr-tenw) 
if (tenw.se.tr )  hln=0. 
c  ceilins-sttic  losses 
3uawe«uawve#araw 
3uaww=uawvw*3raw 
au2c=auc2/295 . 16 
aurn=arn#urfvn 
aurs=ars#urfvs 

3vsum=au2c+3urri+3urs+3uawe+auaww+3c*urfvs 

attmi='=tatic(3uc2Ftern»aurn»  tersraursf  acrurfvsr  tersr  teewr  tewwr 
1  auawe »  auaww  t  svsuni ) 
hlc=au2c# (294 . 16-attmp) 
if ( attmr* . Se . 294 ♦  16 )  hlc=0. 
c  infiltration  losses 
v2=v*2. 30156 
hli*0. 
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if <tat 2t.tr)  hi i=< ♦ 10+3chw*v2+achs* ( tr- ta ) )#vol 
hli=hli#(tr-ta) 
c  window  losses 

wls=u3vs*awins 
if(ta.Se»tr>  wls=0, 
wln=uavn#awinn 
i  f<  ta*s3e*  tr)  wln-O. 
wle=u2ve*awine 
if  <  ta*s*e»tr)  wle^O* 
wlw=usvw*awine 
if  ( ta  ♦  sfe  ♦  tr )  wlw=0. 
win Is- <  wls+wln+wle+wlw ) * ( t r-ta ) 
c  total  loss 

heat  l=hle+hlso+hlw+hln+hlc.+L>ls+hli  twin  Is 

return 

end 
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subroutine  shS ( ansi * rdn*  rduf insul * shsn* mods > 
c  This  subroutine  calculates  the  solar  heat  Sain  throuSh 
c  a  window  as  a  function  of  the  ansle  of  incidence  of 
c  direct  radiation*  the  direct  and  diffuse  components  of 
c  radiation*  and  the  transmissivity  of  the  window, 
c 

dimension  tau2(2) *taud(2) 
dimension  s1p2(2) 
dimension  b2<2> 

data  tau2  < 1 ) » tau2 ( 2 ) /0 . 92  *  0 . 86/ 

data  taud  < 1 ) *  taud  <  2 ) / . 70  * . 62/ 

data  b2(l)*b2<2)/2. 3657*2. 21143/ 

data  s1p2<1) *s1p2<2)/0. 02628 * 0 • 02457/ 

np=insul 

apfc^.SO 

if(np.eo.3)  np=2 
np2=np-l 

if (moda.se. 60i.and.moda.le. 901)  so  to  1 

ans=ansi#57.2957 

if <ans. le*60. )  tauw-tau2<np) 

if  (anS.St  .60.and.anS.  It  .90. )  tauw=sb2(rip)-slp2(np>*ans 
i f ( ans . s t ♦ 90 . )  tauw=0 . 
cans i =cos ( ansi ) 

shsn=< <canSi#tauw*rdn+.5#rdu#t3ud<np) )#(1 ♦ -float ( np2 )# . 10) >*apfc 
So  to  2 

1  shSn=.5#rdu#taud<np) 

2  continue 
return 
end 
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subroutine  fco(v»dir»u3Z> foe) 
c  This  subroutine  calculates  the  surface  film  coeficient 
c  accordinS  to  the  algorithm  given  ba  Kitnura  as  a 
c  function  of  wind  direction  and  speed  relative  to  the 
c  building  surface, 
c 

rwd-was-di r 
rud=abs<  rwd) 

if ( rwd . dt ♦ ISO )  rwd;=360 ♦  -rwd 
if ( rwd .ae .90 . )  do  to  1 
if(v.le.2.)  vc=0.5 
if (v. at. 2)  vcss«25*v 
ao  to  2 

1  continue 
vc=0.3+0.05#v 

2  continue 

foc=67.049*(vc**.605)  +  14.0 

return 

end 


